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ABSTRACT 

We present results on the interstellar medium (ISM) properties of 29 galaxies based on a comparison of 
Spitzer far-infrared and Westerbork Synthesis Radio Telescope radio continuum imag ery. Of these 29 galaxie s, 
18 are close enough to resolve at <1 kpc scales at 70 /jm and 22 cm. We extend the Murphy et alTJ d2006alfbl) 
approach of smoothing infrared images to approximate cosmic-ray (CR) electron spreading and thus largely 
reproduce the appearance of radio images. 

Using a wavelet analysis we decompose each 70 image into one component containing the star-forming 
structures and a second one for the diffuse disk. The components are smoothed separately, and their combi- 
nation compared to a free-free corrected 22 cm radio image; the scale-lengths are then varied to best match 
the radio and smoothed infrared images. We find that late-type spirals having high amounts of ongoing star 
formation benefit most from the two-component method. We also find that the disk component dominates for 
galaxies having low star formation activity, whereas the structure component dominates at high star formation 
activity. 

We propose that this result arises from an age effect rather than from differences in CR electron diffusion 
due to varying ISM parameters. The bulk of the CR electron population in actively star-forming galaxies is 
significantly younger than that in less active galaxies due to recent episodes of enhanced star formation; these 
galaxies are observed within ~10 8 yr since the onset of the most recent star formation episode. The sample 
irregulars have anomalously low best-fit scale-lengths for their surface brightnesses compared to the rest of the 
sample spirals which we attribute to enhanced CR electron escape. 

Subject headings: infrared: galaxies — radio continuum: galaxies — cosmic-rays 



1. INTRODUCTION 

The interstellar medium (ISM) is a complex environment 
comprised of a diverse mix of extremely tenuous matter (by 
terrestrial standards) spanning a range of energetic states. 
Atomic, ionized, and molecular material make up a series of 
gaseous (thermal) phases usually categorized by their tem- 
perature and density; combinations of these components are 
used to define the so-called two- a nd three-phase models of 
the ISM and variations thereof (e.g. lMcKee 1119951 and refer- 
ences therein), however, this is not the entire picture. 

There is an additional phase of the ISM which is often over- 
looked due to the difficulties associated with making direct 
observations of its constituents. This is the relativistic (non- 
thermal) phase which is made up of relativistic charged par- 
ticles known as cosmic rays (CRs) and magnetic fields. The 
CRs within galaxies have an energy density compar able to 
that of the gaseous phases (e.g. lBoulares & Cox I ll990). They 
fill up the entire volume of galaxies and are important sources 
of heating and ionization of the ISM, though characterizing 
thei r propagation remains an ongoing as trophysical problem 
(see[Strong, Moskalenko, & Ptus klnT2007L for a review). 

The relativistic components of the ISM are important dy- 
namically and may play a significant role in the regula- 
tion of star for mation during the formation and evolution 
of galaxies (e.g. iFerriere I l200"lt ICox 1 12005, and references 
therein). Magnetic fields both help to support interstellar 
matter against a galaxy's gravitational potential and con- 
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fine CRs to galaxy disks; thus, magnetic fields and CRs 
take part in the hydrostatic balance and stability of the ISM 
while possib ly determining the properties of gas spiral arms 
(iBeck 1 120071) and even aiding in the triggering of star for- 
matio n jMouschovias. Shu. & Woodward 111974c [Elmegreen 
Il982[). For suf ficiently large CR pressures Parker instabilities 
(Parker"| ["l966l) can create breaches in magnetic disks, allow- 
ing CRs and interstellar material to freely stream into inter- 
galactic space. 

To date, most of our knowledge about the relativistic ISM 
outside of the Galaxy has been obtained indirectly through the 
detection of sync hrotron emiss i on via multi- frequency radio 
observations (e.g. [Duric I 1991; Dahlem, Lisenfeld, & Gollal 
[1991 iLisenfeld et al71[T996l: llrwin et al. I[l99l lB^ck1l2005r) . 
Synchrotron emission arises from CR electron energy losses 
as these particles are accelerated in the magnetic fields of 
galaxies. Although the energy density in CR electrons is 
only ~1% of that for CR nuclei, the similarity between the 
spatial distributions of gamma-ray and synchrotron emission 
within the Galaxy suggests that CR electrons and CR nuclei 
are fairly well mixed on the scales of a few hundred par- 
sees ( e.g. lHaslam et al. | [T982; Bloe men et al. |[l98 6: Webber 
1199 lb . The spatial distribution of a galaxy's synchrotron 
emission is a function of a galaxy's CR electron and magnetic 
field distributions. Thus, radio synchrotron maps provide only 
limited insight on the source distribution of the CR electrons 
as well as the distances the particles may have traveled before 
ending up in their current location of emission. 

Massive stars (>8 M ) are the progenitors of supernovae 
(SNe) whose remnants (S NRs), through the process of dif- 
fusive shock acceleration (Belli [19781: iBlandford & Ostriker I 
1978), appear to be the main acceleration sites of CR 
electrons responsible for a galaxy's observed synchrotron 
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emission. These same young massive stars are often the 
primary sources for dust heating as they emit photons 
which are re-radiated at far-infrared (FIR) wavelengths. 
This shared origin between the FIR and radio emission of 
galaxies is thought to be the foundat ion for the observed 
FIR-radio correlation among (e .g . Ide Jong et al. I 11985 



Helou, Soifer, & Rowan-Robinson 1985; Niklas 



Niklas&Beckl ll997U Yun. Reddy, & Condon 2001) 
Beck & Golla 1 119881; IXu et alT 



within galaxies (e.g 



1997 



and 



1992 



Marsh & Helou 
Hipp elein et al. 
2006). 



1995; Hoernes, Berkhuiisen, & Xu 1998 
2003; Murph y et al. I l2006at iHughes et al 



Coupling the shared origin of a galaxy's FIR and radio 
emission with the fact that the mean free path of dust- 
heating photons (^100 pc) is significantly shorter than the 
expected diffusion le ngth of CR electrons (~l-2 kpc) led 
iBicav & Helou I d 19901) to conjecture that the radio image of 
a galaxy should resemble a smoothed version of its infrared 
image. Consequently, it appears that the close spatial corre- 
lation between the FIR and radio continuum emission within 
galaxies can be used to characterize the propagation history 
of CR electrons. This prescription has been shown to hold 
for galaxies observed at the "super r esolution" (<1') of IRAS 
HIRES data dMarsh & Helou II 1998b and, mo re recently, for 
high re solution (~1 8") Spitzer 70 /jm imaging (iMurphv et al. I 
2006a, hereafter M06a). This phenom enology has been fur- 
ther corroborated on scales >50 pc by Hughes et al. (2006) 
who find synchrotron haloes around individual star-forming 
regions are more extended than FIR-emitting regions within 
th e Large Magellanic Clo ud . 

IMurphv et al. I (l2006bl hereafter, M06b) recently studied 
how the spatial distributions of a galaxy's FIR and radio emis- 
sion vary as a function of the intensity of star formation. They 
concluded that CR electrons are, on average, younger and 
closer to their place of origin within galaxies having higher 
amounts of star formation activity compared with more qui- 
escent galaxies. Using a wavelet-based image decomposi- 
tion, we extend this work by attempting to characterize sepa- 
rately CR electron populations associated with a galaxy's dif- 
fuse disk and its star-forming complexes. We carry out this 
study for a sample of galaxies observed as p art of the Spitzer 
Infrar ed Nearby Galaxies Survey (SINGS; iKennicutt et al. I 
120031) a nd the Westerbork Synthesis Radio Telescope (WSRT- 
SINGS : lBraun et al. 1120071) for which we have the spatial res- 
olution to resolve physical scales < 1 kpc. 

The paper is organized as follows: The galaxy sample is 
defined in §2 while the observations and data reduction tech- 
niques are discussed in §3. In §4 we present and discuss a cor- 
relation analysis between FIR/radio ratios and various other 
physical quantities on sub-kiloparsec scales within galaxies; 
this section leads us to confirm quantitatively that our fa- 
vored phenomenological model is the best description for the 
FIR-radio correlation within galaxies. In §5 we describe our 
wavelet-based, two-component image-smearing analysis; the 
corresponding results are then presented in §6 and their phys- 
ical implications are discussed in §7. We briefly discuss out- 
standing issues and future prospects in §8 and summarize our 
results and conclusions in §9. 

2. GALAXY SAMPLE 

We present FIR and radio continuum imaging for 29 galax- 
ies included in the SINGS (IKennicutt et al. 1 120031) sample. 
The SINGS sample consists of 75 normal galaxies at dis- 
tances <30 Mpc. These objects were chosen to span a range 



of Hubble types (from irregulars to elliptical) as well as ex- 
hibit a large range in star formation rates (SFRs; <0.001 
to >10 M ra yr _1 ), FIR/optical ratios (~10 3 ), and luminos- 
ity (~10 5 ); see Ke nnicutt et al. I ([2003) for a more detailed 
description of the SINGS sample. The WSRT-SINGS sub- 
sample of 29 galaxies was then defined by choosing those 
SINGS galaxies which could be observed with WSRT (i.e. 
North of declination 12°. 5) and have D25 > 5'; these criteria 
ensured that the angular resolution across each target would 
be acceptable to allow for a resolved study of each galaxy 
disk. These 29 galaxies also span a similar ra nge in physi- 
cal par ameters as the parent SINGS sample. See lBraun et al. I 
(2007) for a more detailed description of the WSRT-SINGS 
project. 

In Table[T]we give basic data for each galaxy in the sample. 
Galaxy diameters (D25) and position angles (PA) were taken 
from the Third Reference C atalog of Bright Galaxies (RC3; 
Ide Vaucouleurs. et al. Ill991[). We ca l culate inclinations using 
the method described bv lDale et al. I (1 1997b such that, 
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where a and b are the observed semi-major and semi-minor 
axes and the disks are oblate spheroids with an intrinsic ax- 
ial ratio (b/a)- mt ~ 0.2 for morphological types earlier than 
Sbc and {b/a) m x — 0.13 otherwise. Galaxy distances were 
taken from the literature with preference given to direct mea- 
surements (see Table[T|i. For those galaxies where direct mea- 
surements were not found we use flow-corrected estimates (K. 
Masters 2007, in preparation). Of these 29 galaxies, 18 have 
distances less than 11.5 Mpc; these galaxies are resolved by 
the 70 p.m Spitzer beam at spatial scales <1 kpc. Due to the 
high spatial resolution for which we can examine the FIR- 
radio correlation within these galaxies, they are the focus of 
the present analysis. 

3. OBSERVATIONS AND DATA REDUCTION 

3.1. Spitzer Images 

Observations at 24, 70, and 160 /jm were obtained us- 
ing the Multiband Imaging Photometer for Spitzer (MIPS; 
Rie ke. et al. Il2004l) as part of the SINGS le gacy science pro- 
gram. A detailed descr iption of the SINGS o bservational 
strategy can be found in IKennicutt et al. I (120031) . The MIPS 
data were processed using the MIPS Data Analysis Tool 
(DAT: iGordon et al. Il2005l) and included in the SINGS data 
release 4 (DR4). Additional steps beyond the standard re- 
duction procedure of the MIPS DAT are described in M06a. 
The full width at half maximum (FWHM) of the MIPS 24, 
70, and 160 fjm point spread functions (PSFs) are 5"7, 17", 
and 38", respectively. The mean RMS noise values measured 
for the 24, 70, and 160 p.m maps are ~0.048, 0.447, and 
0.574 MJy sr _1 , respectively. SINGS DR4 has improved the 
final calibration uncertainties to ~4, ~7, and ^12 % at 24, 70 
and 160 /jm, respectively. Some artifacts in the MIPS images 
have also been corrected for in DR4. In the 70 p.m data, very 
high surface brightness regions are known to be affected by 
non-linearities in the calibration. The galaxy most strongly 
affected is NGC 3627 where only ^5 % of the area may be 
miscalibrated. 

3.2. Radio Continuum Images 

Radio continuum imaging at 22 cm was performed using 
the Westerbork Synthesis Radio Telescope (WSRT) as part of 
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TABLE 1 
Basic Galaxy Data 





R.A. 


Decl. 


D 1S 






M opt 


w 10 iw% 


V r 


Dist. 


i 


PA 


Distance 


Galaxy 


(J2000) 


(J2000) 


(arcmin) 


Type 


Nuc. 


(mag) 


(kms" 1 ) 


(kms" 1 ) 


(Mpc) 


o 


O 


References 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


NGC 628 


1 36 41.7 


+15 46 59 


10.5x9.5 


SAc 




-20.9 


74 /175 


657 


7.3 


25 


25 


1 


NGC 925 


2 27 17.0 


+33 34 43 


10.5x5.9 


SABcd 


H II 


-20.6 


224 /267 


553 


9.1 


57 


102 


2 


NGC 2403 


7 36 51.4 


+65 36 09 


21.9x12.3 


SABcd 


H II 


-19.7 


257 /306 


131 


3.2 


57 


127 


2 


Holmb II 


8 19 04.0 


+70 43 09 


7.9x6.3 


Im 




-17.1 


73 /121 


157 


3.4 


37 


15 


3 


NGC 2841 


9 22 02.6 


+50 58 35 


8.1x3.5 


SAb 


Lin/Syl 


-20.7 


611 /664 


638 


14.1 


67 


147 


4 


NGC 2976 


9 47 15.3 


+67 55 00 


5.9x2.7 


SAc 


H II 


-17.6 


.../... 


3 


3.6 


64 


143 


3 


NGC 3031 


9 55 33.2 


+69 03 55 


26.9 x 14.1 


SAab 


Lin 


-21.2 


446/515 


-34 


3.6 


60 


157 


2 


NGC 3184 


10 18 16.9 


+41 25 28 


7.4x6.9 


SABcd 


H II 


-19.0 


142 /396 


592 


11.1 


21 


135 


5 


NGC 3198 


10 19 54.9 


+45 32 59 


8.5x3.3 


SBc 




-20.2 


318 /343 


663 


13.7 


68 


35 


2 


IC 2574 


10 28 21.2 


+68 24 43 


13.2x5.4 


SABm 




-17.7 


123 /134 


57 


4.0 


67 


50 


6 


NGC 3627 


11 20 15.0 


+12 59 30 


9.1x4,2 


SABb 


Sy2 


-20.8 


378 /417 


727 


9.4 


65 


173 


2 


NGC 3938 


11 52 49.5 


+44 07 14 


5.4x4.9 


SAc 




-20.1 


112/265 


809 


13.3 


25 





7 


NGC 4125 


12 08 05.8 


+65 10 27 


5.8x3.2 


E6p 




-21.6 


.../... 


1356 


22.9 


58 


95 


8 


NGC 4236 


12 16 42.1 


+69 27 46 


21.9x7.2 


SBdm 




-18.1 


176 /185 





4.5 


72 


162 


3 


NGC 4254 


12 18 49.5 


+14 24 59 


5.4x4.7 


SAc 




-21.6 


272 /544 


2407 


16.6 


30 





7* 


NGC 4321 


12 22 54.9 


+15 49 21 


7.4x6.3 


SAB be 


Lin 


-22.1 


283 /534 


1571 


14.3 


32 


30 


8 


NGC 4450 


12 28 29.5 


+17 05 06 


5.2x3.9 


SAab 


Lin 


-21.4 


290 /433 


1954 


16.6 


42 


175 


1'' 


rNOC 4JJZ 




. n TO Tl 

+ lz jj Zj 
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Id 
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NGC 4559 


12 35 57.7 


+27 57 36 


10.7x4.4 


SABcd 


HII 


-21.0 


251 /273 


816 


10.3 


67 


150 


7 


NGC 4569 


12 36 49.8 


+13 09 46 


9.5x4.4 


SABab 


Lin/Sy 


-22.0 


360 /397 


-235 


16.6 


65 


23 


7* 


NGC 4631 


12 42 08.0 


+32 32 26 


15.5x2.7 


SBd 




-20.6 


320 /322 


606 


7.7 


83 


86 


10 


NGC 4725 


12 50 26.6 


+25 30 06 


10.7x7.6 


SABab 


Sy2 


-22.0 


410/570 


1206 


11.9 


46 


35 


2 


NGC 4736 


12 50 53.0 


+41 07 14 


11.2x9.1 


SAab 


Lin 


-19.9 


241 /400 


308 


5.0 


37 


105 


8 


NGC 4826 


12 56 43.7 


+21 40 52 


10.0x5.4 


SAab 


Sy2 


-20.3 


311 /363 


408 


5.0 


59 


115 


7 


NGC 5033 


13 13 27.5 


+36 35 38 


10.7x5.0 


SAc 


Sy2 


-20.9 


446 /501 


875 


14.8 


63 


170 


7 


NGC 5055 


13 15 49.2 


+42 01 49 


12.6x7.2 


SAbc 


H II/Lin 


-19.0 


405 /489 


504 


7.8 


56 


105 


7 


NGC 5194 


13 29 52.7 


+47 1 1 43 


11.2x6.9 


SAB be 


H II/Sy2 


-21.4 


195 /244 


463 


7.8 


53 


163 


7 


NGC 6946 


20 34 52.3 


+60 09 14 


11.5x9.8 


SABcd 


HII 


-21.3 


242 /457 


48 


6.8 


32 


69 


11 


NGC 7331 


22 37 04.1 


+34 24 56 


10.5x3.7 


SAb 


Lin 


-21.8 


530/561 


816 


14.5 


71 


171 


2 


References 
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et al. 1 1200111; (5HLeonard et al. 
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NOTE. — Col. (1): ID. Col. (2): The right ascension in the J2000.0 epoch. Col. (3): The declination in the J2000.0 epoch. Col. (4): Major- and minor-axis diameters. 
Col. (5): RC3 type. Col. (6): Nuclear type: H II: H II region; Lin: LINER; Sy: Seyfert (1, 2). Col. (7): Absolute R magnitude, when avai lable; otherwi se from the V 
or B bands. Col. (8): Observed/inclination-corrected 21 cm neutral hydrogen line width at 20% of maximum intensity, in km s~' , taken from Tully 1 1988) or RC3. Col. 
(9): Heliocentric velocity. Col. (10): Distance in Mpc Col. (11): Inclination in degrees. Col. (12): Position Angle in degrees. (13): Distance References 



the WSRT-SINGS survey. A complete description of the ra- 
dio observations an d image processing steps can be found in 
iBraun et al. I d2007l) . To mitigate significant flux loss at low 
spatial frequencies (i.e. the zero-spacing problem), observa- 
tions were taken using a configuration with particularly good 
sampling at short baselines. The WSRT data were CLEANed 
and self-calibrated using an imaging pipeline based on the 
MIRIAD package. The final CLEAN maps were restored 
with 18" FWHM circular Gaussian beams. The intrinsic 
FWHM of the radio beams is approximately 11" East-West by 
1 1/sin (8)" North-South at 1.4 GHz and scales as wavelength, 
where 8 is the source declination. Total intensity calibration 
of the data was performed in the AIPS package and the flux 
density calibration of the radio maps is better than 5 %. In 
the case of NGC 3031, which suffers from a combination of 
low extended surface brightness and confusion problems due 
to the nearby starburst galaxy NGC 3034, a 20 cm map was 
created via a variance weighted average of both 22 and 18 cm 
data in order to obtain high quality data. Accordingly, the flux 
density values of NGC 303 1 where scaled to what is expected 
at 22 cm assuming a mean spectral index of —0.8 to allow 
for proper comparison with the rest of the sample. The mean 
RMS noise in the radio maps is ~-4.4 x 10~ 3 MJy sr _1 . 

3.3. Ancillary Data 

Images in the Ks-hw\A were taken from the Two Mi- 
cron All Sky Survey (2MASS) Large Galaxy Atlas (LGA) 



dJarrett et al 112003b . The typical FWHM of PSF for these im- 
ages is ^2.5". 

Ultraviolet images were obtained by the Galaxy Evolu- 
tion Explorer (GALEX) at 1528 and 2271 A as part of the 
GALE X Atlas of Nearby Galaxies (NGS; iGil de Paz et al. I 
120071) . These galaxies, including nearly all of the SINGS 
sample, were observed with relatively deep integration times 
(~1.5 ksec). With an angular resolution of 4-6", the GALEX 
images are well matched to Sptizer 24 p.m imaging and 
more resolved than the Sptizer 70 and 160 /jm data. For 
a more detailed description of the GALEX observations see 
IGil de Paz et al. I d2007l) . The only galaxy for which we use 
GALEX data that were not obtained as part of the NGS is 
NGC 4725; thes e data were obtained as part of the All-sky 
Imaging Survey dMartin et al. 1 120051 AIS;) and have much 
shorter effective exposure times (^0.1 ksec). Galaxies in 
the present sample which have not yet been observed by 
GALEX, or for which the data are not publicly available, in- 
clude NGC 3 1 84 and NGC 6946. 

3.4. Image Registration and Resolution Matching 

We match the resolution of the MIPS and radio images us- 
ing Gaussian PSFs rather than the MIPS PSFs, which suf- 
fer from significant power in their side-lobes. After cropping 
each set of galaxy images to a common field of view and re- 
gridding them to a common pixel scale, we divide the Fourier 
Transform (FT) of the final MIPS images at each band by the 
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FIG. 1 . — Each galaxy's 70 /.im and 22 cm images are displayed in columns 1 and 2, respectively. These images are shown using a logarithmic stretch ranging 
from the 1 -G RMS level of the background to the maximum surface brightness of the galaxy. A normalized 570 map and the residual images between the 22 cm 
and 70 /jm images, smoothed with the b est-fi t single smearing kernel and best-fit disk and structure smearing kernels are given in columns 3, 4 and 5 respectively. 
The residual map definition is given in 95.41 The stretch of each of these 3 maps runs from -0.75 to 0.75 dex; red and blue colors correspond to radio and infrared 
excesses in the residuals, respectively, while green corresponds to residuals ~0. Regions removed for the residual calculations (e.g. background radio sources) 
appear as dark, circular holes in a few of the residual maps. 



FT of a model of the corresponding MIPS PSF. Next, this quo- 
tient is multiplied by the FT of a Gaussian PSF matching that 
of the CLEAN beam used to restore the corresponding radio 
images. The final product was then taken back into the im- 
age plane where we checked to ensure that flux was conserved 
throughout the resolution matching procedure. The resolution 
matched 70 fjm and 22 cm maps for each galaxy are presented 
in first and second columns of Figure Q] 

The 2MASS A^-band images were, likewise, registered, 
convolved, and re-gridded to match the MIPS 70 /jm data. 
The same was done for the GALEX far-UV (1528 A; FUV) 
and near-UV (2271 A; NUV) images after each was first cor- 
rected for Milky Way ex t inctio n using the Galactic color ex- 
cesses of ISchlegel et al. I d!998l) . The FUV and NUV images 
were then added to create a single GALEX "total UV" image 



for each galaxy which we have used in the present analysis. 

Once the PSF matching of the registered 70 and 22 cm 
maps was completed, qjo maps were constructed where, 



q 10 = log 



/v(70 A< m)[Jy] 
S v (22 cm)[Jy] 



(2) 



3.5. Global Flux Densities 

Global flux densities at 22 cm and at each MIPS band were 
measured by summing the flux density at each wavelength 
within elliptical apertures. Any identifiable flux contributions 
from background sources were masked out in the measure- 
ments. In the case of NGC 5194, emission from its compan- 
ion galaxy was also excluded in the global flux density mea- 
surements. These MIPS and radio flux densities agree with 
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FIG. 1. — Continued 



those presented in lDale et al~l (120071) and lBraun et al. I (12007b 
(within errors). We therefore choose to adopt their values for 
consistency purposes. The only notable exceptions are the 
22 cm flux densities of NGC 3031 and NGC 5194. Our ra- 
dio flux density value for N GC 303 1 w a s mea sured using the 
global aperture defined by iDale et aT71 d2007l) and our aver- 
aged 20 cm map (see 83.2| >. We then scaled this flux density 
to what is expected at 22 cm assuming a mean spectral index 
of —0.8 to allo w for proper compa rison with the rest of the 
sample. Unlike iBraun et al. I d2007l) . we mask out flux con- 
tributions from NGC 5195 when measuring the 22 cm flux 
density of NGC 5194. The monochromatic 22 cm, 24, 70, 
and 160 /jm flux densities for each sample galaxy are listed in 
Tabled 

Using a weighted combination of the MIPS flux densities 
we compute estimates of the total-infrared (TI R; 3 — 1 100 fim) 
flux of each galaxy according to Equation 4 of Da le & Helou I 
(2002). We also compute an estimate of the IRAS based far- 
infrared (FIR; 42 — 122 /jm) fractions using the same spectral 



energy distribution (SED) models of IDale & Helou I (12002b 
such that 

(3) 

^tir £6 V/v(160//m)y 

where [^^1,^2^3] = [0.5158,0.1619,-0.3158,-0.1418] 
and -0.65 < log / v (70//m)// v (160 jan) < 0.54. 

The TIR fluxes were then used to calculate global ^tir ra- 
tios for each galaxy where, 

F TlR \ , / S v (22cm) 



q-m = log 



The more commonly used q 

iHelou. Soifer. & Rowan- Robinson (1985), which 
denote as, qpm is also calculated such that 

F F1R \ , ( S v (22cm) 



3.75 x 10 12 Wm- 2 



commonlv 



log 



W m 2 Hz 



used 



(4) 
of 

we 



<?FIR = log 



3.75 x 10 12 Wm- 2 



log 



Wm~ 2 Hz 1 



(5) 



6 



MURPHY ET AL. 



i 4 



E«>dnjJ R«>d 



NQC 4736 



* # # 



NQC 4326 



NQC 5055 



'^flfe^ j^pP^ 

p 

NQC51S4 

# ^ ^ 



NQC 



Fig. 1. 



Continued 



The FIR and TIR fluxes, along with the associated q^in and 
<?tir ratios, are given in Table [2] 

3.6. Radiation Field Energy Densities 

We compute the radiation field energy density (£/ ra d) of each 
galaxy using its TIR surface brightness since this parameter 
is sensitive to the diffusion of CR electrons. This was done 
by first repeating the above image registration and resolution 
matching at each MIPS waveband using a 45" Gaussian beam 
(i.e. slightly larger than the FWHM of the native 160 /jm 
beam) creating TIR surface brightness maps for each galaxy 
using a weighted combination of the 3 M IPS band according 
to Equation 4 of iDale & Helou I ([2002) on a pixel by pixel 
basis. The TIR luminosity (Xtir) was then summed within 



,-2 <,„.-! 



elliptical apertures fit to a galaxy's 1.4 x 10 W m sr 
isophotal radius; this radius corresponds to the maximum 3-o 
RMS value among the entire sample's TIR maps and was used 
to calculate the isophotal TIR-diameters Dtir for each galaxy. 
We note that these luminosities will be less than those corre- 



sponding to the TIR fluxes given in Table [2] as those fluxes 
were measured within much larger apertures, eff ectively go- 
ing int o the noise of each MIPS map, as given by Dale et alTI 
(2007). Using the deprojected area of each elliptical aperture, 
Ajir = 7i(Dtir/2) 2 , we calculate TIR surface brightnesses 
(Etir = ^tir/Atir) along with estimates of t/ ra d for radia- 
tion emitted near the surface of a semi-transparent body such 
that 



t/rad ; 



2Jt Ltir 

'bol Z 7T A 

c 2A T irc 




(6) 



where /boi is a galaxy's bolometric surface brightness, c is the 
speed of light, and all quantities are given in cgs units. The 
parenthetical term in Equation|6]provides a correction fo r non- 
absorb ed UV emission that was empirically derived by Bell 
(2003) using archived FIR and UV data for a sample of more 
than 200 galaxies. 
Since f/ ra d is likely proportional to the disk-averaged star 
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TABLE 2 

Global Flux Densities and Derived Parameters 





i v (2z cm) 


4-' /*") A . .m \ 

/ v (24/jm) 


/ v (70/jm) 


/v(160 fjm) 




"fir 




na 






Galaxy 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(10~ 13 WmT 2 ) 


(10~ 13 Wm~ 2 ) 


<?FIR 


<?TIR 


NOC 628 


0.200 


3.19 


34.78 


126.2 


17.36 


49.39 


2.36 


2.82 


wr,r Q9^ 

J *i v 1 v . i 


090 


95 


14 40 


43 33 


6.67 


17 47 


2 30 


2 71 


NGC 2403 


0.360 


5.84 


86.36 


245.6 


39.64 


101.5 


2.47 


2.88 




005 


20 


3 67 


4 46 


1.36 


2 73 


2 82 


3 12 


Nnr 9X41 

L \ V. J V. > Z.OtJ 


100 


91 


10 22 


62 29 


6.53 


20 80 


2 24 


2 74 


l \ v. j v . z.y i u 


068 


1 37 


20 43 


52 56 


9.18 


22 61 


2.56 


2 95 




590" 


5 09 


85 18 


360 


41 83 


128 5 


2 28 


2 77 


NGC 1 1 84 


0.080 


1.43 


15.76 


70.48 


8.10 


25.69 


2.43 


2.93 


Nnr ^iqs 


049 


1 06 


10 27 


39.00 


5 24 


15 24 


2 46 


2.92 


IC 2574 


0.014 


0.28 


5.55 


11.75 


2.32 


5.33 


2.65 


3.01 


Nnr 1^97 

IN vH JOZ / 


n sno 


7 42 


Q9 f\1, 


LJXJ.L 


49 24 


102 8 


2 35 


2 74 


Nnr iq^r 

L \ V. J V. J7JO 


080 


1 09 


14.25 


51 98 


6.97 


19 87 


2 37 


2 82 


Nnr 4.17s 


002 


079 


111 


1 77 


0.46 


96 


2 81 


3 13 


Nnr 47^6 


026 


55 


8.27 


20 43 


3.67 


8 92 


2 58 


2 96 


NGC 4254 


0.510 


4.20 


50.29 


142.9 


23.69 


60.65 


2.09 


2.50 


NGC 4321 


0.310 


3.34 


40.59 


139.6 


19.79 


54.95 


2.23 


2.67 


NGC 4450 


0.013 


0.21 


3.42 


16.94 


1.82 


5.79 


2.57 


3.08 


NGC 4552 


0.093 


0.094 


0.54 


1.42 


0.29 


0.72 


0.92 


1.31 


NGC 4559 


0.110 


1.12 


16.89 


54.15 


7.93 


21.35 


2.28 


2.71 


NGC 4569 


0.170 


1.44 


12.37 


41.21 


6.29 


17.23 


1.99 


2.43 


NGC 4631 


1.290 


8.15 


130.2 


289.5 


56.22 


131.5 


2.06 


2.43 


NGC 4725 


0.100 


0.86 


8.85 


59.91 


6.17 


19.66 


2.22 


2.72 


NGC 4736 


0.320 


5.65 


93.93 


177.4 


38.94 


86.47 


2.51 


2.86 


NGC 4826 


0.110 


2.72 


55.16 


98.82 


22.06 


48.28 


2.73 


3.07 


NGC 5033 


0.240 


1.97 


28.81 


91.07 


13.54 


36.21 


2.18 


2.60 


NGC 5055 


0.450 


5.73 


72.57 


302.3 


36.43 


111.1 


2.33 


2.82 


NGC 5194 


1.240* 


12.67 


147.1 


494.8 


71.84 


197.5 


2.19 


2.63 


NGC 6946 


1.700 


20.37 


207.2 


502.8 


97.10 


234.1 


2.18 


2.56 


NGC 7331 


0.590 


4.36 


74.97 


189.5 


33.01 


80.79 


2.17 


2.56 



References. — Radio flux densities taken from Braun et al. 1 120071) : MIPS flux densities taken from Dale et all 120071) . 
NOTE. — a : This flux density, measured using a 20 cm map, was scaled to the expected flux density at 22 cm assuming a spectral index 
of —0.8. b : Corrected for flux contributions from companion galaxy NGC 5195. 



formation activity within the normal star-forming galaxies 
considered here, we also express t/ ra d, for illustration pur- 
poses , as a star fo rmation rate surface density, Esfr, again 
using iBelil d2003l) . The TIR diameters along with values of 
Stir and t/ ra d, expressed in units of L fi kpc~ 2 and erg cm~ 3 , 
respectively, are given in Table[3] 

4. CORRELATION ANALYSIS OF q 70 

M06a showed that the FIR/radio ratio is not constant within 
galaxies, exhibiting strong variation patterns in residual maps 
and a dispersion as small as the dispersion in the FIR/radio 
ratios among galaxies. Some of the variation in the FIR/radio 
ratios within galaxies was found to be related to both infrared 
surface brightness and radius. Using our larger sample, we 
confirm these results here as well as expand the explored pa- 
rameter space. We search for the dominant physical parameter 
related to the variations in the FIR/radio ratio within galaxies 
by plotting qjQ against a number of physical quantities and 
calculating the residual dispersion in each fit. This was done 
for galaxies sufficiently nearby that the resolution at 70 is 
better than 1 kpc. 

4.1. q-jQ- Correlation Analysis within Galaxies 

We examine the relations between q-/Q and: a) 24 /jm sur- 
face brightness: since it ha s been found to trace local massive 
star formation activity (e.g. lCalzetti et al. I200ll2007h : b) Ri- 
band surface brightness: tracer of stellar mass; c) UV surface 
brightness: extinction dependent tracer of local star forma- 
tion; d) galactocentric radius; e) 24/70 /jm surface brightness 
ratio: local heating of dust; f) 24 /jm/UV surface brightness 
ratio: obscured star-formation and permeability of ISM. Each 



TABLE 3 

Radiation Field Energy Densities and Star 
Formation Rate Surface Densities 



Dtir log Stir log U mA (TIR+UV) 
Galaxy (kpc) (Z, kpc -2 ) (erg cmT 3 ) 



NGC 628 


20.42 


7.37 


-12.67 


NGC 925 


24.42 


6.94 


-13.06 


NGC 2403 


12.79 


7.38 


-12.60 


Holmb II 


4.76 


6.57 


-12.91 


NGC 2841 


30.69 


7.21 


-12.86 


NGC 2976 


5.98 


7.48 


-12.38 


NGC 3031 


20.12 


7.19 


-12.82 


NGC 3184 


21.33 


7.42 


-12.63 


NGC 3198 


27.07 


7.16 


-12.89 


IC 2574 


10.68 


6.34 


-13.32 


NGC 3627 


23.10 


7.80 


-12.30 


NGC 3938 


20.85 


7.48 


-12.57 


NGC 4125 


13.79 


6.88 


-13.00 


NGC 4236 


15.50 


6.33 


-13.44 


NGC 4254 


31.47 


7.80 


-12.31 


NGC 4321 


30.07 


7.67 


-12.43 


NGC 4450 


20.17 


7.11 


-12.89 


NGC 4552 


4.51 


6.70 


-12.82 


NGC 4559 


24.52 


7.14 


-12.89 


NGC 4569 


21.50 


7.59 


-12.48 


NGC 4631 


27.57 


7.58 


-12.51 


NGC 4725 


30.59 


7.01 


-13.03 


NGC 4736 


13.76 


7.61 


-12.42 


NGC 4826 


8.14 


7.77 


-12.20 


NGC 5033 


40.23 


7.25 


-12.84 


NGC 5055 


26.45 


7.56 


-12.52 


NGC 5194 


30.46 


7.68 


-12.42 


NGC 6946 


27.81 


7.73 


-12.38 


NGC 7331 


44.81 


7.47 


-12.64 
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TABLE 4 
Residual Dispersion Results 



Galaxy 


< <?70 > 


070 


a c [/ v (24^m)] 


Oc [Radius] 


a c [/ v (24^m)// v (70pm)] 


a c [/ v (&)] 


o c [/ v (UV)] 


C c [/ v (24pm)// v (UV)] 


NGC 628 


2.26 


0.258 


0.166 


0.189 


0.258 


0.183 


0.182 


0.221 


NGC 925 


2.24 


0.205 


0.170 


0.200 


0.190 


0.211 


0.177 


0.196 


NGC 2403 


2.27 


0.253 


0.133 


0.171 


0.248 


0.165 


0.115 


0.238 


NGC 2976 


2.36 


0.344 


0.153 


0.248 


0.344 


0.207 


0.204 


0.316 




2.22 


0.290 


0.212 


U.248 


0.236 


U.23U 


n net 
0.256 


0.259 


NGC 3184 


2.22 


0.226 


0.173 


0.204 


0.202 


0.195 






NGC 3627 


2.07 


0.219 


0.179 


0.187 


0.203 


0.187 


0.188 


0.209 


NGC 4236 


2.31 


0.454 


0.395 


0.452 


0.447 


0.515 


0.411 


0.446 


NGC 4559 


2.23 


0.232 


0.159 


0.206 


0.224 


0.199 


0.170 


0.213 


NGC 4631 


1.76 


0.245 


0.204 


0.226 


0.237 


0.183 


0.221 


0.242 


NGC 4725 


2.24 


0.266 


0.155 


0.247 


0.237 


0.216 


0.241 


0.230 


NGC 4736 


2.14 


0.234 


0.206 


0.208 


0.193 


0.210 


0.208 


0.234 


NGC 4826 


2.45 


0.211 


0.183 


0.179 


0.167 


0.173 


0.169 


0.209 


NGC 5055 


2.06 


0.244 


0.196 


0.214 


0.168 


0.213 


0.206 


0.212 


NGC 5194 


1.92 


0.234 


0.172 


0.210 


0.233 


0.176 


0.175 


0.232 


NGC 6946 


1.93 


0.256 


0.179 


0.235 


0.244 


0.228 






Averages 


2.17 


0.261 


0.190 


0.226 


0.239 


0.218 


0.209 


0.247 



of the UV, Ks, 24 /jm, 70 /jm, and 22 cm images were re- 
gridded with 1 kpc 2 pixels. Pixels in the re-gridded images 
of each galaxy were removed if they were either below the 3- 
a level or associated with emission from non-related objects, 
such as background galaxies. Galaxies having <10 such res- 
olution elements remaining were excluded in the analysis (i.e. 
Holmb II and IC 2574). 

An example of the qjo scatter-plots for each of the six tested 
parameters are given for NGC 628 in Figure [2] Each rela- 
tion is fit using an ordinary least squares (OLS) regression of 
the form qjQ = mlog (X) +b, where X is one of the six pa- 
rameters cited above, m is the slope of the fit, and b is the 
y-intercept of the fit. To quantitatively determine which of 
these six parameters best correlates with the FIR/radio ratio, 
we compute the residual dispersion about the regression line, 
denoted as Oc [X] . Average qjo ratios and the associated dis- 
persions are given in Table H] along with the residual disper- 
sions for each parameter and galaxy; the means for the entire 
sample are also presented. 

4.2. Results of q-/o Correlation Analysis 

We find that the residual dispersion around the trend 
of qjQ versus 24 fim surface brightness is, on average, 
^0.07 dex smaller than the dispersion in qjo', this is the 
largest improvement among all the parameters inspected. 
The average values of Gc[Iv(Ks)], Gc[/v(UV)], <Jc [Radius], 
a c [/ v (24 ^m)// v (70 /an)], and a c [/ v (24 /mi)/7 v (UV)], are 
-0.03, 0.02, 0.04, 0.05, and 0.05 dex larger, respectively, 
than the average value of Oc[A>(24 fan)] suggesting that the 
FIR/radio ratio is most strongly correlated on 24 /jm sur- 
face brightness. A galaxy's 24 /jm surface brightness appears 
to be a good proxy for extinction corrected Paa and there- 
fore the amount of ongoing massive star formation activity 
within a galaxy (e.g. ICalzetti et al. Il200"5h : however, we note 
that a combination of 24 /jm and H a emission, to account 
for unobscured star formation, has bee n found to be a better 
proxy than 24 fjm em ission alone (i.e. ICalzetti et al. 1 120071: 
iKennicutt et al. Il2007l) . The fact that the residual dispersion 
between qjQ and UV surface brightness is nearly as small as 
the residual dispersion around the trend of qjQ and 24 /jm sur- 
face brightness suggests that, at least at kpc scales, most of 
the star formation within our sample is not highly ob s cured . 
This result is consistent with that of iPrescott et al. I (2007) 
who find that highly-obscured star forming regions make up 



only ~3% of a sample of ^1800 star-forming regions within 
SINGS galaxies on scales down to ^500 pc. 

Comparing the parameters of radius and Ks surface bright- 
ness, we find that the residual dispersions are similar. This 
result is not too surprising as the Ks surface brightness sam- 
ples the older stellar population of a galaxy which, generally, 
has a smooth spatial distribution as a function of radius. We 
also find that the residual dispersion around the trend of qjo 
versus / v (24 /jm)// v (70 /jm) is only slightly smaller than the 
dispersion in qjQ, and significantly larger than the residual dis- 
persion around the trend of qjo versus 7 V (24 fjm) suggesting 
qjo is not very sensitive to temperature changes. 

4.3. Local FIR/Radio Ratio and Star Formation 

The above results demonstrate that, at least on kpc scales, 
the local FIR-radio correlation is strongly coupled to star- 
formation activity. However, from an empirical standpoint, 
identifying physical parameters to minimize the residual dis- 
persion in qn) is not where the greatest improvement has been 
found. While we find an average decrease in the residual dis- 
persion between qjQ and 24 fjm surface brightness of ^20%, 
M06b has shown that a phenomenological smearing model 
improves the correlation between the FIR and radio emission 
distributions of galaxies by a factor of ^1.8. This signif- 
icantly larger improvement using the image-smearing tech- 
nique likely arises because this method more precisely charac- 
terizes the time evolution of a galaxy's FIR and radio emission 
distributions, which depend critically on its star formation his- 
tory. The image-smearing procedure minimizes the residuals 
between a galaxy's infrared and radio images after smooth- 
ing the infrared map by a parameterized kernel assumed to 
contain the diffusion and decay characteristics of a galaxy's 
CR electrons. Since CR electrons spread over time, the size 
of each smoothing kernel is a measure of diffusion and decay 
which captures the time dependence mentioned above. Other 
parameters (e.g. density, magnetic field strength, etc. . . ) can 
affect the permeability of the ISM to CR electrons and impede 
their propagation; such parameters will vary with CR electron 
location, from the vicinity of star-forming complexes to the 
more diffuse regions of the ISM. As the image-smearing pic- 
ture appears to provide the best quantitative description for the 
FIR-radio correlation, we will now introduce a more sophisti- 
cated approach to take into account the differences among CR 
electron populations within galaxy disks. 
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FIG. 2. — Plots comparing the dependence of ^70 on 24 ,um surface brightness, Ks surface brightness, GALEX UV surface brightness (FUV + NUV), 
galactocentric radius, [7 V (24 /im)// v (70 fim)] flux density ratio, and [7 V (24 J um)// v (UV)] flux density ratio are plotted in panels (a), (b), (c), (d), (e), and (f), 
respectively, for NGC 628. Each data point corresponds to a measurement from 1 kpc 2 pixels. Ordinary least squares (OLS) fits, which were used to measure the 
residual dispersion, are overplotted in each panel. 



5. TWO-COMPONENT IMAGE-SMEARING ANALYSIS 

Using a p henomenologica l imag e-smearing model, first 
presented by lBicay & Heloul (I1990I) . M06b studied how the 
spatial distributions of FIR and radio emission varied as a 
function of star-formation intensity within 12 spiral galaxy 
di sks. Their result s were then related to the physical model 
of iHelou & Bicavl ([T993). Although this method was found 
to improve the spatial correlation between the distributions 
of infrared and radio emission by reducing the residuals in 
the ratio maps, a single kernel failed to work perfectly for 
an entire galaxy disk as evidenced by the recognizable struc- 
tures (i.e. spiral-arms) left in the residual maps. This is 
not too surprising as the freshly injected CR electron popu- 
lations associated with star-forming regions should be differ- 
ent in mean age and energy than those populations associated 
with a galaxy's diffuse disk; such CR electrons have likely 
lost a significant amount of their initial energy and/or have 
undergone re-acceleration by passing interstellar sho cks (e.g. 
IStrong & Moskalenko |[l99l iBerezhko et al. 1120031) . Conse- 
quently, it appears that a more realistic description should re- 
quire a multi-scale analysis of the galaxy images allowing for 
the separation of structures at various spatial frequencies. We 
now present such an analysis using a wavelet-based image- 
decomposition technique and compare these results to those 
using a single-component method. 

The simplest way to formulate a phenomenological model 
with two distinct CR electron populations is to create repre- 
sentative source functions and allow them to be characterized 
by different scale-lengths. We therefore decompose the in- 
frared images (used as source proxies) into two components; 
one representing sites of ongoing star formation and the other 
consisting of diffuse emission from dust in the disk heated by 
older stars or from afar by young stars. While a purely spatial 
definition will not yield a perfect census of massive star for- 
mation sites, it has the advantage of simple assumptions and 
data manipulation. 



5.1. Choice of Separation Scale 



The goal of the decomposition is to separate each 70 
image into two components: one containing the smallest to 
largest physical structures associated with massive star for- 
mation and another to include a galaxy's diffuse FIR emission 
distribution. The typical distance for which >90% of Lyman 
continuum photons are ex pected to be absorbe d by interstellar 
gas and dust is ^500 pc (Dove & Shulll [19941) : this value is 
in excellent agreement with the observed m ean-distances of 
diffus e ionized gas surrounding HII regions (Ferguson et al. 
1996). One kpc is also a moderate upper-limit to the sizes of 
giant HII regions and molecular cloud complexes, as well as 
widths of spiral arms. Furthermore, the results presented by 
M06b hinted that the best-fit global scale-lengths seemed to 
be spread around a median value of ~1 kpc, but did not cluster 
tightly around it. For these reasons we choose a disk-structure 
separation scale of 1 kpc. While this scale may not be perfect 
for all galaxies, it is physically motivated and allows for a uni- 
form comparison among the nearest 18 WSRT-SINGS sample 
galaxies. 

5.2. Wavelet Based Image Decomposition 

The wavelet transform provides the ability to describe struc- 
tures within an image as a function of their characteristic spa- 
tial scale. This property of the wavelet transform has led to the 
develop ment of the wavelet based Multi-scale Vision Mode l 
(MVM; iBiiaoui & Ruil [1991 iRue. F. & Biiaoui. A.I [1991 . 
a procedure useful for ide ntifying morphological features in 
astron omical images (e.g. IStarck et al. 1 120001 ; lAdami et al. I 
120051) . Wavelet analyses have also been used in to study 
the FIR-radio correlation within individual systems such 
as NG C 6946 (Fric ketal. Il200l. the LM C iHughes et al.1 
(2006), and M 33 (Tabat abaei et al. 112007b : in each of these 
studies a cross-correlation analysis was performed on the 
wavelet power spectra for images acquired at various wave- 
lengths. 

For the purpose of this study, we need only make use of the 
wavelet transform to separate structures in each image as a 
function of their characteristic spatial frequency. We apply the 
wavelet transform to our images using the a trous algorithm 
dHoldschneider et al. lll989tlBiiaoui 11199 ll) . The a trous algo- 
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rithm produces a set of wavelet coefficients w p , at each scale 
p, having the same number of pixels as the original image, 
c„. The exact physical scale corresponding to each wavelet 
plane, as indexed by p, will depend on the distance to the ob- 
ject and the pixel scale of the image. It then follows that a 
pixel j of the observed image can be expressed as the sum all 
the wavelet coefficients at this position, plus a final smooth 
plane containing the remaining large scale features, cp, such 
that, 

p 



CoJ = cp J + 



E w p 

p=i 



(7) 



An example of the wavelet transform via the a trous algorithm 
is given in Figure [3] for NGC 6946. The complete a trous 
algorithm, as well as detailed discussions on the applicability 
of wavelet transforms to astr onomical data sets, can be found 
in lStarck & Murtaghl d2002l) . 

We decompose each observed infrared image, /, into two 
sub-images: (1) a structure image containing features with 
spatial scale smaller than 1 kpc; (2) a disk image containing 
all structures with characteristic spatial scales larger than or 
equal to 1 kpc, largely constituting a galaxy's diffuse disk. 
The a trous algorithm is dyadic, transforming the image at 
discrete scales which are powers of 2 and indexed here by 
p. We therefore re-grid each image appropriately such that 
2 p ~ l pixels, where p > 3, corresponds to 1 kpc at a given 
galaxy's distance. This re-gridding scheme ensures that the 
structure image is sampled at the Nyquist frequency and pro- 
vides uniformity in the decompositions for all galaxies. The 
final structure and disk images, 7 str and 1^, respectively, are 
then created by summing the appropriate scales such that 



p(<\ kpc) 

E 

p=i 



and 



^d S k=^p+ E w p> 

p{>\ kpc) 



(8) 



(9) 



where Ip is the final smooth plane containing all remaining 
large scale features and / = 7 str + /dsk- For illustration, we dis- 
play in Figure [4] the structure and disk images at 70 pm for 
three galaxies (NGC 628, NGC 5055, and NGC 2976) chosen 
to span a large fraction of our sample's range in infrared sur- 
face brightness. Each panel is displayed with the same stretch 
to allow comparison between the morphology and amount of 
power residing in the disk and structure components as a func- 
tion of infrared surface brightness. 

5.3. Correcting for Free-Free Emission 

The image-smearing analysis is designed to compare the 
spatial distributions of a galaxy's FIR and non-thermal radio 
emission. Because the observed radiation at 22 cm is a com- 
bination of non-thermal (synchrotron) and thermal (free-free) 
emission, we estimate and subtract the thermal fraction of ra- 
dio emission from the observed 22 cm maps. This is done 
by estimating the free-free emission using a scaled version of 
24 pm maps as described in M06b. 

The scaling factor relies on the empirical correlation found 
to exist between 24 pm and extinction corrected Paa lumi - 
nosities within NGC 5194 (M51a) by ICalzetti et al. I (120051) . 
While it ha s been demonstrated that this correlation is not uni - 
versal (e.g. Perez-Gonzalez et al. ll2006t[Calzetti et al. 120071) . 



TABLE 5 
Radio Thermal Fraction 
Comparison 



f l GHz 

Jt 



fl GHz 

Jt 



Galaxy Niklas, Klein, & Wielebinsky (1997) This Paper 



NGC 3031 


0.05 


0.06 


NGC 3627 


0.09 


0.09 


NGC 4254 


0.05 


0.05 


NGC 4321 


0.05 


0.06 


NGC 4569 


0.07 


0.07 


NGC 4631 


0.03 


0.04 


NGC 4736 


0.18 


0.10 


NGC 5033 


0.11 


0.05 


NGC 5055 


0.09 


0.08 


NGC 5194 


0.05 


0.06 


NGC 6946 


0.06 


0.07 


NGC 7331 


0.06 


0.06 



it has proven to be more than sufficient as a first order es- 
timate of the thermal radio emission (M06b). Having esti- 
mated the Paa line emission from the observed 24 pm emis- 
sion, we can then determin e the corresponding ionizing pho- 
ton rate (Osterbrock 119891) and expected free-free emission 
at 1.4 GHz <|RubjnJll968|) such that, 



- 7.93 x 10 



0.45 



-0.1 



3 I * \ I \ - f M24pm) \ 

iq 4 k) vghz/ V J y / 

(10) 

where we have assumed an average H II region temperature 
ofT = 10 4 K. 

To check the reliability of this method for estimating 
the radi o thermal fractions, we compare our results to 
those of iNiklas, Klein, & Wielebinsky I i 1997b for 12 galax- 
ies whi ch appear in both samples. The thermal frac- 
tions of INiklas. Klein. & Wielebinsky I d 19971) were derived 
by the completely independent method of radio spectral 
index fitting. We list the 1 . GH z thermal fractions of 
INiklas. Klein. & Wielebinsky! dl997l) . along with our esti- 
mated thermal fractions, after scaling our 1 .4 GHz flux densi- 
ties to what is expected at 1 .0 GHz assuming a mean spectral 
index of —0.8 for the non-thermal component, in Table|5] We 
find that, except for the cases of NGC 4736 and NGC 5033, 
our estima ted thermal fractions are very simil ar to those re- 
ported by lNiklas. Klein. & Wielebinsky I d!997l) . We note that 
NGC 5033 hosts a Seyfert-2 type nucleus which could com- 
plicate this method of free-free estimation using the 24 pm 
imaging. 

5.4. Phenomenological Image-Smearing Model 

The basic procedure used here is similar to that presented 
in M06a,b. We calculate the residuals between the free-free 
corrected radio and observed infrared images after convolv- 
ing the infrared maps by a parameterized kernel K(r). The 
new element is that we now smear 7 str and /dsk indepen- 
dently, then add the two smeared images and compare the 
sum to the radio image. The smoothing kernel is a function 
of a two-dimensional position vector r, having a magnitude 
r — (x 2 +y 2 y/ 2 , where x and y are the right ascension and 
declination offsets on the sky. 

Due to the large range in inclination among our sample 
galaxies, we compare results for exponential smoothing ker- 
nels oriented in either the plane of the galaxy disk or isotropi- 
cally. An exponential kernel was chosen because it was found 
to work as well as, or better than, Gaussian kernels by M06a. 
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FIG. 3. — The wavelet transform of NGC 6946 at 70 fjm using the a trous algorithm (see i|5.2l . In this instance the image was re-gridded such that the each 
wavelet plane w p= i-s corresponds to 2 P x 62.5 pc and the final smoothed plane, cp contains informatio n for the remaining larger spatial scales. For the actual 
creation of the structure and disk images used in the present analysis, the re-gridding scheme described in j)5.2l was used which ensures that the structure image is 
sampled at the Nyquist frequency. Each plane is displayed using a logarithmic stretch ranging from the l-o~ RMS level of the background to the maximum value 
of the wavelet coefficients at that plane. The observed 70 /jm image is recovered exactly by summing these six planes. 



This is consistent with "leaky box" type models where CR 
electron escape is expected to occur on tim escales less than 
or comparable to the diffusion timescales (Bic av & Heloul 
1990). The kernel takes the form, K(r) = e~ r l r °, where r Q 
is its e-folding length / modified by geometric factors related 
to the position angle 9 and inclination ; of a galaxy such that 



/ cos i 



[1 - (xsine + ycos9) 2 sin 2 «7r 2 ] 1 /2' 



(11) 



The position angle of the tilt axis of the galactic disk is 
measured East of North and / = defines a face-on projection. 
Kernels are isotropic when r Q = /. This situation is the same 
as if a galaxy's inclination and position angle were equal to 
zero. 

Let R(r) denote the observed radio image and /(r) = 
Atr(r) +Aisk( r ) denote the observed infrared image, which 
has been separated into component images containing struc- 
tures larger or smaller than ~1 kpc (see ^5. 2b . The two- 
dimensional residual function, defined by (j), is calculated be- 
tween the radio and smeared infrared images after the infrared 
disk and structure images are first smoothed by independent 
kernels and then summed. This residual surface is given by 



<KG, Wdsk) 



(12) 



where Q — L^/( r )/L^;( r ) is used as a normalization factor 
[i.e. log Q = qx(global)], 



J(4tr,'dsk) =4tr(4tr) +^dsk(^dsk) 



(13) 



represents the infrared image after smearing 7 str and /^sk with 
independent kernels having scale length / str and i<jsk, respec- 
tively, and the subscript j indexes each pixel. Our calculation 
of the residuals is slightly different than what has been used in 



M06a,b; we now fix the normalization factor, Q, rather than 
let it vary with the kernel size. This modification is necessary 
since we are now using different smoothing kernels for each 
component. The quantity (]) was minimized to determine the 
best-fit smearing kernel for the structure and disk images of 
each galaxy in our sample. The normalization by the squared 
sum of the radio flux density allows for proper comparison of 
our galaxies which vary in intrinsic surface brightness. We 
plot log (j) as a function of / str and i<jsk in Figure[5] The special 
case of a single kernel, as studied by M06a,b, is reproduced 
by the diagonal cut of of the residual surface plotted in Figure 
|5](i.e. when / str = /dsk)- These residual curves are presented in 
Figure [6] and binned by the amount of star formation activity 
within each galaxy. 

The estimation of the residuals is carried out after first re- 
moving pixels that are not detected at the 3-a level in either 
the radio or maximally-smeared infrared maps, and editing 
out contaminating background radio sources. Additional edit- 
ing was necessary for NGC 3031, NGC 3627, NGC 4725, 
NGC 4736, and NGC 4826. Each of these galaxies host 
an AGN. Because the phenomenology we are trying to test 
is strictly associated with the physics of star formation pro- 
cesses, we try to remove any contamination an AGN might 
introduce. This is done by first subtracting a scaled PSF from 
the nuclear region of the 70 /jm images before smoothing to 
remove any excess emission associated with the AGN. We 
then excise a 1 kpc region around the nucleus before calculat- 
ing the residuals between the smeared-infrared and radio im- 
ages. Since the AGN within SINGS galaxies are all very low 
power, a 1 kpc region should be adequate to remove the in- 
fluence of an AGN. Additional editing was also necessary for 
NGC 5194; its companion galaxy (NGC 5195) was removed 
before calculating the residuals. 
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FIG. 4. — The structure, disk, and total (observed) 70 /jm images of NGC 628, NGC 5055, and NGC 2976. These galaxies were chosen to span our sample's 
range in infrared surface brightness and are ordered, left to right, by increasing radiation field energy density. We have used a common logarithmic scaling for all 
six panels running from the 3-G level in the noisiest map to the maximum surface brightness level among all the maps. The structure image contains all featu res 
significant at spatial frequencies less than 1 kpc while the disk image contains all features significant at spatial scales greater than or equal to 1 kpc. (See j|5.2l for 
details.) 

images as measured by the l-o RMS noise of each map; the 
uncertainty on the best-fit scale-length is then estimated as the 
range in scale-length corresponding to the range from min((|)) 
to min((j))+unc((|)), spanning along the diagonal and structure 
and disk axes of the residual surface maps. 

While (j) and <t> characterize the residual behavior as sin- 
gle quantities, we also constructed residual images for each 
galaxy using the best-fit smearing kernels to inspect the spa- 
tial variations of the residuals using this new, two-component 
method. The residual images, defined by, 

Residual image = log (<2 _1 f(r)) - log (R(r)), (16) 

are plotted in column 5 of Figure Q] for isotropic kernels hav- 
ing an exponential profile as this kernel type was associated 
with the largest values of <f> (see $63). For a comparison, the 
residual maps associated with the best-fit global scale-lengths 
are given in the fourth column of Figure Q] 



We use the quantity, 

<S> = logf ♦(g.O.O) V (14) 
Vmin[(j)(2,/str,/dsk)]/ ' 

which is the depth of the minimum value in the residual sur- 
face, as a measure of how much the correlation is improved 
by smoothing the infrared structure and disk images. In the 
special case when / str = Zdsk, we have 

= log ( ■ jJn ^ i m ) ' (15) 

\mm[^(Q,l str = l dsk )\J 

which is the same <t> studied by M06a,b (i.e. the minima found 
for the diagonal residual curves presented in Figure |6j. We 
denote the best-fit scale-length which maximizes Ogic* as / g i b 
and refer to this quantity as the best-fit global scale-length for 
the galaxy. Uncertainties in (]) are estimated by numerically 
propagating the uncertainties in the input 70 um and 22 cm 
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FIG. 5. — The residual surfaces for each galaxy which are ordered by increasing U m( i from left to right, top to bottom. The abscissa and ordinate designate the 
scale-length of the kernel used to smooth the structure and disk component images, respectively. The minimum value on the surface is identified with an asterisk. 
A diagonal cut (i.e. where / str = 1^) through each surface plot corresponds to the one dimensional residual curves which are plotted in figure[6] 
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FIG. 6. — Res iduals between the radio maps and smeared 70 /im images (as 
defined in 35.41 as a function of smearing scale-length. The results plotted 
here use isotropic kernels and 22 c m ma ps which have been corrected for the 
presence of free-free emission (see j|5.3t . Galaxies having low disk-averaged 
star f ormation rates, defined by log U m i < —12.5 = log £sfr < —2.3 (see 
j|3.6l are plotted in panel (a) while galaxies with high star formation activity, 
and larger values of l/ rac | (£sfr)> are plotted in panel (b). 



An additional free parameter of the relative amplitude of the 
disk and structure components was considered for a subset of 
the sample spanning the full range in surface brightness. This 
was done to test whether the proposed two-component analy- 
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FIG. 7. — Radiation field energy density binned by morphological type. 
The plotting symbol sizes are scaled as a function of the improvement found 
using the two-component image-smearing procedure compared to using the 
single-component method (i.e. — O^i,). This value ranges from 0.00 dex 
for NGC 4736 to 0.24 dex for NGC 463 1 . 

sis alone is sufficient at capturing the signatures of CR elec- 
trons associated with a galaxy's structures and disk or whether 
the relative efficiencies between structure and disk processes 
in producing synchrotron emission is significant. This could 
be especially important in the cases for which the power in the 
structure and disk components are comparable. While slight 
quantitative improvements (i.e. on the few percent level) are 
seen, they were not dramatic enough to warrant extra free pa- 
rameters. 

6. TWO-COMPONENT MODELING RESULTS 

While we have shown that the FIR/radio ratios measured 
within galaxies are most sensitive to star formation activity in 
§4, we note that the largest quantitative improvement to the 
dispersion in the FIR/radio ratio within galaxies is achieved 
via our image-smearing analysis (i.e. nearly a factor of 2 com- 
pared to a factor of ~ 1 .2 using a residual dispersion analysis); 
this procedure describes the spreading and decay of CR elec- 
trons in galaxy disks. In fact, by repeating the correlation 
analysis described in §4 with the best-fit smoothed infrared 
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maps we find that the gradients seen between q-j® and 24 /jm 
surface brightness and radius are flattened, on average, by fac- 
tors of 5 and 4, respectively. Hence, our phenomenological 
model clearly captures the behavior observed in §4; infrared 
emission is more strongly peaked around star-forming regions 
since CR electrons diffuse away from bright structures into 
diffuse regions of the ISM. 

In the following section we present the results from our 
two-component image-smearing analysis and look for the 
cases in which significant gains are obtained over the single- 
component method. Since this phenomenology is valid for 
non-thermal radio continuum emission, we focus our discus- 
sion on the results using radio images that have been corrected 
for the presence of free-free emission as described in jj5.3l and 
note that this correction only slightly increases (i.e. 0.1 kpc 
on average) the measured best-fit scale-lengths. 

6.1. Comparison with Single-Component Method 

We find a range in the behavior of the residual surface maps 
among the sample galaxies in Figure [5] For nearly all galax- 
ies we observe an initially steep (negative) gradient for in- 
creasing structure scale-lengths as the minimum in log (j) is 
reached. Even though some galaxies may not exhibit a well 
defined minimum along the structure axis (e.g. NGC 628), 
the steep initial gradient clearly demonstrates that the struc- 
ture component is being detected by our two-component anal- 
ysis. Similarly, many galaxies show a well defined minimum 
along the disk axis (e.g. NGC 4559) while others do not (e.g. 
NGC 4631) suggesting that this component is evident as well. 
This diverse behavior among the residual maps demonstrates 
that our two-component analysis does yield qualitative insight 
into the two major CR electron populations under discussion. 

To assess quantitatively the gains of using the two- 
component approach over image-smearing with a single ker- 
nel, we look at the differences between <t> and Ogiob which 
we give in Table [7] We find a large range in improvement 
spanning from 0.00 (NGC 4736) to 0.24 dex (NGC 4631). 
While galaxies showing the largest improvements generally 
have high infrared surface brightnesses, galaxies exhibiting 
very little improvement are well distributed as a function of 
infrared surface brightness. We also assess whether improve- 
ment correlates with morphology. 

In Figure [7] we plot each galaxy's radiation field energy 
density as a function morphological type; the plotting symbol 
sizes are proportional to 4> — <b % \ h- Irregular and early-type 
spirals show essentially no improvement by using the two- 
component approach. Significant improvement appears lim- 
ited to spiral galaxies of type Sb and later with high surface 
brightnesses. 

The negligible improvement found for irregulars (4> — 
«t>gi b <0.005 dex) can be explained by their lack of a diffuse 
emission component at both FIR and radio wavelengths. This 
is clearly seen by inspecting the 70 pm and 22 cm images of 
these galaxies in Figure Q] As a consequence of this, these 
galaxies have been excluded in the statistical comparisons of 
best-fit scale-lengths presented in §6.3 and §6.5; these galax- 
ies will be revisited in jj7.3l of the discussion. 

As for the early-type spirals, small (<t> — Ogiob ^0.02 dex) 
improvements are probably expected since their FIR and ra- 
dio emission are dominated by an inner-disk; if the cen- 
tral concentration of emission largely dominates over that 
of the disk, this morphology will be pretty well described 
by a single-component compared to a galaxy disk contain- 
ing many organized star-forming structures. Consequently, 



we find that the two-component approach has mainly im- 
proved our phenomenological description of the FIR-radio 
correlation for late-type galaxies exhibiting clear spiral struc- 
ture and significant star-formation throughout their disks; this 
corresponds to roughly half the sample spirals, for which 
< 4> 4> g iob >~0.1 dex. 

A discrepancy is found for NGC 2976, which is a mem- 
ber of the M 8 1 group. This discrepancy may be due to the 
fact that the galaxy's optical morphology (which is classified 
as an Sc peculiar) is strikingly different than its morphology 
in the FIR. In the FIR its disk emission arises from a very 
small disk (Dtir ^6 kpc) and is dominated by two large star- 
forming complexes on either side of a relatively weak nu- 
cleus; these features are almost unnoticeable in the optical. In 
fact, detailed studies of NGC 2976 carried out in the optical 
(Bron kalla. Notni. & Mutter] 1 19921) and in HI (IStil & Israel I 
12001 " have shown that the galaxy has dwarf like character- 
istics and possibly triggered star formation due to encoun- 
ters between other members in the central M8 1 group. Fur- 
thermore, the optical luminosity, and presumably mass, of 
NGC 2976 is much smaller than that for the sample spirals, 
but similar to the sample irregulars. Such findings are con- 
sistent with our result; NGC 2976 behaves more like the ir- 
regulars in the sample for which we find little improvement 
is obtained by using the two-component analysis compared to 
the single-component method. 

6.2. Morphological Differences in Residual Images 

The residual images using the single-component and two- 
component analyses are displayed in the fourth and fifth 
columns of Figure [T] respectively. As already stated in 36.11 
we find that the decrease in the residuals between the FIR and 
radio maps for the two-component analysis with respect to 
the single-component analysis is significantly greater for late- 
type spirals having a relatively large amount of ongoing star 
formation. 

The residual images using the single-component method 
display a general trend such that galaxies with low star for- 
mation activity exhibit radio excesses associated with star- 
forming regions while galaxies with higher star formation 
activity behave in an opposite manner; these galaxies dis- 
play infrared excesses associated with star-forming regions in 
their residual images. This is the result of galaxies with low 
star-formation activity generally needing large scale-lengths 
to best match their infrared and radio spatial distributions 
leading to over-smoothed star-forming regions relative to the 
more actively star-forming galaxies. Now, by using the two- 
component technique and smoothing small spatial scale fea- 
tures (i.e. individual star-forming complexes) with smooth- 
ing functions independent of those used for the large spatial 
frequency components, this trend has been largely been sup- 
pressed for a number of galaxies (i.e. NGC 4631, NGC 5194, 
NGC 6946). 

A comparison of the residual images for our sample's near- 
est grand-design spiral (NGC 6946) in Figure Q] best illus- 
trates this result; the rather dramatic differences in the resid- 
ual behavior for star-forming spiral arms and quiescent inter- 
arm regions, which are evident when using a single kernel, 
are largely suppressed by our new treatment. Some orga- 
nized structure still remains in our two-component residual 
image which is not surprising since this method only looks 
to characterize the extremes of a galaxy's CR electron pop- 
ulation. The remaining structures might also suggest that 
our symmetric smoothing functions have limited applicabil- 
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FIG. 8. — In the bottom portion of each panel we plot radial cuts across the 70 fjm and 22 cm maps of NGC 925, NGC 4631, and NGC 5194. The 22 cm profile 
has been scaled by a factor of 100 to allow for a better comparison against the 70 /.im profile. In the top portion of each panel we plot the ratio between the radio 
and the observed infrared profiles along with the ratio between the radio the infrared profiles after smoothing by the best-fit kernels. Each cut, displayed in the 
thumbnail image of the galaxy, has a width equal to the FWHM of the 70 fjm PSF. 
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Effects of Inclination on Kernel Shape for Sample Spirals 
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ity; diffusion preferentially occurring along field lines, which 
trace arms and become tangled in star-forming regions, will 
lead to asy mmetries asso ciated with spiral arm structure (e.g. 
NGC 6946: lBeck1l2007h . 

Since it is difficult to determine amplitude information from 
the residual images in Figure [1] we show radial cuts across 
the 70 fjm and 22 cm maps, along with the ratios between 
the radio and both the observed and smoothed 70 /jm im- 
ages, in Figure [8] This is done for three galaxies (NGC 925, 
NGC 4631, and NGC 5194) chosen to span our sample's mor- 
phology and inclination. These plots clearly illustrate that the 
amplitude of the variations in the FIR/radio ratios across each 
galaxy disk is significantly reduced by our image-smearing 
procedure. In the case of the grand design spiral NGC 5194, 
we also find that the arms are broader in the radio compared to 
the FIR; this observation is consistent with CR electrons being 
accelerated in star formation sites within spiral arms and dif- 
fusing larger distances than associated dust heating photons. 

6.3. Inclination Effects 

We now determine how inclination affects the choice of 
smoothing kernel. In Table|6]we give the median values of <I> 
and Ogiob for isotropic or in-disk kernels after separating the 
sample into low (i < 60°) and high (i > 60°) inclination bins. 
As expected, the difference for face-on galaxies appears al- 
most negligible; the median SO = med(<I>i S0 — Ogai) for these 
galaxies is ~0.02 dex while SOgiob is also ^0.02 dex. Con- 
versely, if a galaxy's inclination is greater than ^60° then the 
orientation of the kernel seems to become important. In this 
case, isotropic kernels are strongly favored such that SO is 
~0.17 dex while 80 g i b is ~0.13 dex. 

The regular component of a galaxy's magnetic field is 
spread most densely throughout its thin disk; the radial dif- 
fusion of CR electrons should then preferentially occur along 
field lines while vertical (out-of-plane) diffusion should re- 
quire CR electrons to undergo an increased amount of cross- 
field diffusion. Such a scenario has been verified empiri- 
cally. Indirect estimates of diffusion coefficients for the ver- 



tical propagation of CR electrons in a galaxy's thin disk have 
been found to b e an order of magnitude smaller tha n those for 
radial diffusion (iDahlem. Lisenfeld. & Goilall 19951) . 

However, in galaxies with active star formation, ordered 
magnetic fields can be ruptured allowing CR electrons to 
quickly escape their disks and form synchrotron haloes (e.g. 
iHummel et al. lfT988h . This decrease in vertical confinement 
will lead to a diffusion behavior which is more isotropic in 
appearance. A clear example of this scenario is seen in the 
nearly edge-on galaxy, NGC 4631. NGC 4631 has a large 
radio halo that extends ~5 kpc beyond the vert ical extent 
of the FIR disk dDahlem. Lisenfeld. & Golia1[l995T) . Our re- 
sults therefore suggest that the highly inclined sample galax- 
ies each possess, at least to some degree, synchrotron haloes 
in which the diffusion of CR electrons occurs on similar time- 
scales as those in the disk. 

6.4. Best-fit Scale-lengths 

In Figure [5] we plot the residual surfaces associated with 
the two -component image-smearing procedure as described 
in {35.41 The best-fit disk, structur e, an d global scale-lengths, 
along with estimated errors (see {35.41 i. are given in Table [7] 
We find that the best-fit scale-lengths for the structure im- 
ages are on average ~ 0.5 kpc while the best-fit scale-lengths 
for the disk images are, on average, ~1.3 kpc; or, in other 
words, the best-fit disk scale-length is, on average, ^2.5 times 
larger than the best-fit structure scale-length. The best-fit 
global scale-lengths, using the single-component approach, 
are found to be ~ 1.0 kpc, on average. 

To ensure that the exact separation scale used in decompo- 
sition of the 70 /jm images does not affect the best-fit disk and 
structure scale-lengths we repeated our analysis using sepa- 
ration scales ranging from 0.5 (for galaxies resolved at this 
scale) to 2.0 kpc for 4 galaxies of varying size, morphology, 
and infrared surface brightness (i.e. NGC 628, NGC 2976, 
NGC 2403, NGC 6946). We find that the dispersion in best- 
fit disk and structure scale-lengths is ~0.15 and 0.10 kpc, re- 
spectively. Our separation method therefore seems quite ro- 
bust and the best-fit disk and structure scale-lengths appear to 
be fairly insensitive to the exact separation scale. 

In Figure |9h we plot the best-fit disk and structure scale- 
lengths versus radiation field energy density, which we also 
express as a star formation surface density. Also included in 
this plot are the best-fit scale-lengths found using the single- 
component method. We perform ordinary least squares fits to 
the best-fit disk, structure, and global scale-lengths with f/ rat j 
as the independent variable. The least square fits for the best- 
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fit disk and structure scale-lengths in Figure^ have very sim- 
ilar slopes (-0.43 and -0.47, respectively). Furthermore, the 
best-fit global scale-lengths appear to be nearer to the best- 
fit disk scale-lengths for galaxies having low values of t/ ra d 
and nearer to the best-fit structure scale-lengths for galaxies 
having high values of {/ ra d- This result is consistent with an- 
other observation, namely that the best-fit global scale-lengths 
display a general trend of decreasing with increasing ratio be- 
tween the power in the 70 structure to disk components. 
These results suggest that the difference between the FIR and 
radio emission spatial distributions is dependent on whether 
the disk or structure is the dominant emission component of a 
galaxy. In §7 we propose that the disk will generally be com- 
posed of a population of CR electrons which are significantly 
older than those which have been recently injected into the 
ISM and still reside close to their parent star-forming regions. 

We also note that four galaxies (NGC 925, NGC 4725, 
NGC 4559, and NGC 3031) appear to be off of the gen- 
eral trend between best-fit global scale-length and star for- 
mation rate surface density first identified by M06b. While 
the best-fit structure scale-lengths of these galaxies appear 
consistent with the rest of the sample, the best-fit disk and 
global scale-lengths appear to be lower than what one would 
expect compared to the trend for the higher surface bright- 
ness galaxies. Of the 15 galaxies we plot in Figure [9^, these 
four galaxies also have the four lowest surface brightnesses. 
This turnover at logf/ ra d ~ —12.7 may be due to a signal-to- 
noise effect; specifically, we may not be detecting the most 
diffuse emission in the extended parts of these galaxy disks 
which is throwing off our estimations for the best-fit scale- 
lengths. For example, by truncating the areas over which we 
calculate the residuals, the turnover is shifted towards values 
of log£/ ra( i > —12.7. Conversely, simulations of much deeper 
observations to detect a very low surface brightness disk lead 
to longer scale lengths, suggesting again that the turnover may 
be due to the limited sensitivity of our data. On the other 
hand, the location of these galaxies in Figure^ may not be a 
signal-to-noise artifact, but due to a different physical make- 
up of these systems. Because it is difficult to settle this issue 
without deeper observations, we include a discussion on the 
physical implications of this scenario in 37.31 

7. DISCUSSION 

We find that the FIR and non-thermal radio morphologies 
are more similar to each other for galaxies having higher radi- 
ation field energy densities compared to galaxies with lower 
radiation field energy densities. Following the interpretation 
of M06b, our results indicate that CR electrons are, on aver- 
age, closer to their place of origin in galaxies having higher 
star formation activity. 

Since the diffusion scale-length of CR electrons depends 
only on their age and ability to diffuse through the ISM of 
galaxies, there are four possible explanations of the results. 
The CR electrons within galaxies having high star formation 
activity may: (1) have relatively short lifetimes due to a high 
energy loss rate; (2) diffuse at a slower rate due to the ISM 
having a high density and magnetic field strength, resulting in 
a shorter mean free path; (3) escape into intergalactic space 
at a higher rate due to an increased occurrence of ruptured 
magnetic field lines (this explanation implies systematically 
higher global FIR/radio ratios); or (4) have been accelerated 
recently and be relatively young. The first three of these ex- 
planations are applicable in the case of steady-state star for- 
mation (i.e. the variation of ISM parameters alone will lead 



TABLE 7 
Two-Component Fitting Results 





'glob 






^glob 




3> — ^lob 


Galaxy 


(kpc) 


(kpc) 


(kpc) 


(dex) 


(dex) 


(dex) 


NGC 0628 


n i +0.73 
Z,A -6.§1 


U - J -0.18 


<"j t+0.60 
• -042 


0.499 


0.561 


0.063 


NGC 0925 


0.8jjHg 


£+0.37 
U -°-0.20 


i 1+0.57 
1 - 1 -0.42 


0.479 


0.501 


0.022 


NGC 2403 


i - u -0.21 


„ _-i_n ifi 
u ■ J -o.io 


^ -0.32 


0.664 


0.762 


0.098 


Holmb II 


u,A -0.04 


n . +0 10 
A 1 ^ 
u - 1 -0.04 


~ ™4-1 44 


0.043 


0.043 


0.000 


NGC 2976 


u -^-o.oi 


u -^-0.03 




0.864 


0.867 


0.003 


NGC 3031 


i q+0.73 
LS -0.58 


a ^+0.59 
U.J_ 022 


i a+0.66 
1 -7-0.44 


0.464 


0.503 


0.039 


NGC 3184 


. ~4-n 39 
L -*-0A5 


a Ai-U-O/ 
uy -0.26 


z - u -0.35 


0.741 


0.764 


0.023 


IC 2574 


a n+0.17 
u - z -0.03 


A O+O.ll 

u -'-0.13 


A A+0.60 

u - u -NaN 


0.161 


0.167 


0.006 


JNGC 362/ 


a c+0.01 


A ,1+0.01 

°- 4 -0.03 


A Q+0.13 

US -0.07 


0.541 


0.629 


0.088 


NGC 4236 


„ -4-fi Of, 

l\ /■) -|-u.zo 

u - z -0.16 


4-0 I 1 

ul -0.12 


.4-3 11 


0.009 


0.017 


0.008 


NGC 4559 


, ,+0.03 
'•-Q.QS 


a fi +0.06 
uo -0.11 


, ,,+0.06 
1 ■*-<). 14 


0.655 


0.713 


0.058 


NGC 4631 


(, - 6 -0.09 


4+ " 1 
u -^-o.oo 


, ,+0.05 
1 - 1 -0.06 


0.768 


1.009 


0.241 


NGC 4725 


, 1+0.84 
1 " 5 -0.47 


a n+3.94 
uy -0.41 


1 6 +0M 
1,o -0.71 


0.570 


0.585 


0.015 


NGC 4736 


u -°-o.oi 


6 +0 03 
u -°-0.04 


ofi +0.11 
u -°-0.07 


0.729 


0.729 


0.000 


NGC 4826 


a 9+0.01 
u - z -0.04 


a 9+0.01 
u - z -0.04 


5 +0 - 23 
u - J -0.25 


0.117 


0.138 


0.021 


NGC 5055 


U - S -0.09 


5+ 01 
u - -0.00 


, ,+0.02 
1 - -0.03 


0.727 


0.805 


0.078 


NGC 5194 


A Cj+0.10 

u " 3 -0.04 


A -+0.01 

u -^-0.05 


, ,+0.13 
1 - 1 -0.14 


0.343 


0.417 


0.073 


NGC 6946 


6+0.14 
u -°-o.io 


a 4 +0.02 

u -^-o.oo 


1 6+ 010 
10 -0.11 


0.690 


0.860 


0.171 


Averages" 


o.y 


0.5 


1.3 


0.5% 


0.656 


0.066 



NOTE. — °: The irregular galaxies (Holmb II, IC 2574, and NGC 4236) were ex- 
cluded when calculating the averages. 



to shorter scale-lengths). The fourth explanation, however, 
requires a recent episode of enhanced star formation to inject 
fresh CR electrons into a galaxy's ISM. We will try to distin- 
guish among these four scenarios in the context of the results 
for the two-component smearing analysis. 

7.1. CR-Electron Cooling Timescales and Diffusion 
Properties in Normal Galaxies 

We now describe the physical processes associated with the 
propagation and decay of CR electrons in the ISM. We de- 
rive CR electron cooling timescales and introduce a diffusion 
equation to derive propagation distances and physical scaling 
relations for comparison with our observations. These com- 
parisons will allow us to discriminate between the four sce- 
narios above. 

As CR electrons propagate through the ISM of galax- 
ies they lose their energy due to a number physical pro- 
cesses including synchrotron, inverse-Compton (IC) scat- 
tering, bremsstrahlung, ionization, and adiabatic expansion 
losses. In normal galaxies synchrotron and IC scattering pro- 
cesses are likely the most significant energy l oss terms for CR 
electrons associated with 1 GHz emission (ICondonlll992h : 
the other terms listed will become non-negligible, however, 
for galaxies hosting extreme episodes of star-formation like 
starbursting ultra-lum inous infrared galaxies (ULIRGs) (e.g. 
Thompson et al. 2006). We now derive the CR electron cool- 
ing timescales associated with the two dominant processes. 

Let us assume that CR electrons propagating with a pitch 
angle a in a magnetic field of strength B have isotropically 



distributed velocities such that < sin a >= I leading to 
B± w 0.82 B. According to synchrotron theory, a CR elec- 
tron having energy E will emit most of its energy at a critical 
frequency v c where 

2 
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Since the energy loss of CR electrons by synchrotron radiation 
goes as dE /dt°c UbE 2 , we can use EquationfTTlto express the 
synchrotron cooling timescale, x syn = E j\dE /dt\ syn , for CR 
electrons as 

_ ~ 1/2 (JL) 1/2 ( Ub 

GHz/ \fjGJ ViO^^ergcm- 3 

(18) 

where Ub = B 2 /{%%) is the magnetic field energy density. 
Naturally, we can rewrite the synchrotron cooling timescale 

as 



^M~5.7xl0 7/ V 
yr 



L syn 

yr 



1.4 x 10 



(—) 
VGHz/ 



1/2 ( B 
vi uG 



-3/2 



(19) 



Similarly, since the energy loss of CR electrons due to IC scat- 
tering goes as dE/dt U^E 2 , we can again use EquationfTTI 
to write the IC cooling timescale, Tic = E/\dE/dt\ic, as 



Tic 
yr 



<5.7x 10 



7 f— ) 
VGHz/ 



-1/2 



B 



1/2 



rad 



10 12 erg cm 3 
(20) 

where t/ rac i is the radiation field energy density. Photons with 
a frequency v p will significantly contribute to IC losses if be- 
low the Klein-Nishina limit; for GeV electrons considered 
here, the upper photon energy limit lies in the X-ray regime 
(i.e. /zv p < 200 eV), but the bulk of losses arise from interac- 
tions with IR/optical photons, which dominate C/ la( j. The ef- 
fective cooling timescale for CR electrons due to synchrotron 
and IC losses is 

Tic 1 , ( 21 ) 



-l -l 

T cool — % syn ' 



which, by combining Equations [T8l and l20l we can express as 



Tcool 

yr 



<5.7x 10 



VGHz/ 



-1/2 



B 



U2 



u B + u, 



10 12 erg cm 3 
(22) 

In simple diffusion models, the propagation of CR elec- 
trons is usually characterized by an em pirical, energy- 
depen dent diffusion coefficient, De (e.g. iGinzburg et al. I 
1980). The value of De has been found to be around 4 — 
6 x 10 28 cm 2 s _1 for <GeV CRs by fitting diffusion models 
with direct measurements of CR nuclei (i.e. secondary-to- 
primary ratios li ke Boron-to-Carb o n) within the Solar Neigh- 
borhood (e.g. [J ones et alTl l2001t iMoskalenko et al. 1 120021: 
iMaurin. Taillet. & Donato 1 120021) . While this empirically 
measured value is for CR-nuclei within the Milky Way, it 
has been found to be consistent with inferred diffusion co- 
efficients for CR electrons both radially (~ 10 29 cm 2 s ) 
and vertically (~ 10 28 cm 2 s -1 ) in the thin disks of galaxies 
(i.e. N GC 891 and NGC 463 1 : iDahlem. Lisenfeld. & Golla I 
1 1995b . We note that a value of ~ 10 29 cm 2 s _1 is also found 
for hydrodynamic simulations of bubble/super-bubble in- 
duce d galaxy outflows in s t udies of galaxy cluster abundances 
(e.g. iRebusco et al. Il200l iRoediger et al. 112007b . Since the 
mean model-derived diffusion coefficient, using direct mea- 
surement of CR nuclei, is similar to the mean value of in- 
fer red radial and vertical CR e l ectron diffusion coefficients 
bv lDahlem. Lisenfeld . & Golla] d 19951) . we simply use 



D E 



5 x 10 28 , 

5xl0 28(^)l/2, 



E < 1 GeV 
E > 1 GeV. 



(23) 



Now, neglecting escape and using a simple random- walk 
equation, we assume CR electrons will diffuse a distance 



'cool = (De t C oo1 )'^ 2 before losing all of their energy to syn- 
chrotron and IC losses. By combining Equations 1171 and l23l 
we can express De as a function of B for a fixed v c such that, 
for CR electrons having energies > 1 GeV, 



'cool 

kpc 



7x 10~ 



(if) " 2 ( 



GHz 



B 



-1/8 



(24) 



7.1.1. Order- of -Magnitude Estimates 



Using the above equations, we derive simple, order-of- 
magnitude estimates to determine whether diffusion and cool- 
ing of CR electrons in a steady-state star formation model 
are able to account for our observations. In the following 
sections, more careful analyses are presented. Taking the 
mean value of £/ ra d for those galaxies plotted in Figure |9j3 
(i.e. 3.7 x 10~ 13 erg cm -3 ), and assuming £/ ra( j = Ub, which 
has been shown to be a reasonable assumption for a large 
'sample of spiral galaxies (Lisenfe ld. Volk. & Xul[l996l) . we 
find from Equation [22] that the average cooling time for a 
1.4 GHz emitting CR electron is ~1.1 x 10 8 yr. Inserting this 
value into Equation [24] we measure a diffusion scale-length 
of ^6.8 kpc; this value is clearly off of the scale shown in 
Figure [5J). On the other hand, if we inst ead assume a fixed, 
typical magnetic field strength of 9 /jG dNiklas 111 9951) . and 
that Ub = C/ r ad = 3.2 x 10~ 12 erg cm~ 3 , the average cooling 
time for a 1 .4 GHz emitting CR electron is ~2.2 x 10 7 yr with 
a diffusion scale-length of ^2.6 kpc. While this value for t/ ra d 
is much higher than what we infer from the average TIR sur- 
face brightness of the sample, it must apply near bright star- 
forming structures, whose TIR surface brightnesses are much 
_ 1 greater. Even so, this diffusion scale-length is much larger 
than any value we find for the best-fit structure scale-lengths. 
From these simple order of magnitude estimates it appears 
that particle fading due to cooling by Inverse Compton and 
synchrotron processes alone cannot explain the structural dif- 
ferences between the FIR and radio maps. 

While CR electron escape may help to reconcile this dis- 
crepancy between the observed scale-lengths and these order 
of magnitude diffusion scale-length calculations, we might 
then expect to find systematic variations in global FIR/radio 
ratios with best-fit scale-lengths. No such trend is found, 
suggesting that escape is not a dominant variable among the 
galaxies included in Figure [9}}. The role of escape, however, 
may be important for the sample irregulars which is discussed 
in 37.31 For the sample spirals, the best explanation seems to 
be differences in the CR electron population ages. 

7.1.2. Comparison with Scaling Relations 

We now introduce a more detailed model to see if we can re- 
produce the observed trends in Figure [9}} by varying relevant 
ISM parameters using the above mentioned scaling relations. 

Again setting Ub = U a d leads 5 toS« t/ r ' ad 2 and T C ooi x B^l 2 . 
Since De Z m fp and /mf p °= « 5 , where Z m f p is the mean-free- 
path traveled by CR electrons, n is the ISM density, and 8 
depends on the specifics of the scattering processes within a 
system, we can estimate the relative importance of B and n 

5 We note that this relation between B and U m d is slightly steeper than the 

B t* C/ r ' a j 3 scaling reported by[Niklas & Beck] U997T) . If we were to instead 
use this scaling, the final relation in Equation l26l would be less steep and even 
more discrepant with the fit to our observations in Figure|9J>. 
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FIG. 9. — In panel (a) we plot the best-fit disk, str uctu re, and global scale-lengths for galaxies in our sample which are resolved at scales less than 1 kpc 
against the radiation field energy dens ities, (C/ la( i, see 33.61 - We exclude the galaxies Holmb II, IC 2574, and NGC 4236 which have morphologies not well fit 
by our phenomenological model (see j|6.2> . Least square fits for the best-fit disk, structure, and global scale-lengths are p lotte d as dot-dash, dashed, and dotted 
lines, respectively. In panel (b) we have excluded the galaxies NGC 925, NGC 4725, NGC 4559, and NGC 3031 (see §6A\ . Along with the fit to the global 
scale-lengths (dotted line) we also plot with the expected diffusion scale-lengths due to inverse Compton (IC) losses in a fixed mag netic field (long-dashed line) 
and synchrotron + IC losses with an energy-dependent diffusion coefficient Dg for the steepest possible index (solid line; see j|7.1t . 
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are realistic index ranges 
ence 



on the distance CR electrons emitting at a characteristic fre 
quency travel. We introduce the scaling relations 

BocnP, D E ocn 5 B~ 1/4 

where ± <P< |and- l < 8< -5 
dHelo u & Bicavl [T993l) . thereby leading to the dependi 

7 „ .,1/4(8/13-7/4) 

'cool x U rad . (26) 

While the range o f p has been derived empirically (e.g. 
iNiklas & Beck 1 1 19971) . and is consistent with the assump- 
tion of equipartition between turbulent kinetic energy of gas 
clouds a nd magnetic energy, as well as with dyn amo mod- 
els (e.g. iRuzmaikin. Sokoloff. & Shukurovl 119881) . the val- 
ues of 8 are determined by the dominant CR electron scat- 
tering processes. A value of —1 describes how diffusion 
would occur for hard scattering off particles (i.e. the hard 
sphere approximation). On the other hand, 8 = — J best de- 
scribes soft scattering off of structures; specifically, the case 
in which scattering is dominated by localized B perturbations 
whose spatial density i s proportional to the density of the gas 
dHelou & Bicav |[l993l) . In reality, both scattering processes 
are likely active which will correspond to an intermediate 
value of 8; this value will likely vary within and among galax- 
ies. 

Excluding those galaxies which may suffer from low 
signal-to-noise effects, we first take the simplest case and con- 
sider whether the observed trends in Figure^ may arise from 
an increase of IC losses alone within galaxies that are more 
actively forming stars. Fixing B, we plot the expected relation 

between Z coo i and t/ ra d [i.e. / coo i (< U B > +(7 ra d) -1 / 2 ] as a 
dashed line in Figure |9j?. The fit to the global, structure, and 
disk data have slopes that are ~7.7, 2.7, and 5.1 times steeper 
than the slope of this line. It is probably more physical for 
each parameter, not just C/ ra d, to scale with the amount of star 
formation activity within a galaxy; we will now investigate 
such a case. 

Using the derived scaling relation given in Equation|26] and 
taking p = i and 8 = — 1, we plot the expected trend between 
'cool and t/ rat j in Figure |9]3 as a solid line. These extreme 
choices of p and 8 correspond to the steepest possible slope 
and are unlikely to best describe most galaxies. To more eas- 
ily compare this model with the fits to the disk, structure, and 
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FIG. 10. — The slopes of the fits to global, structure, and disk scale-lengths 
as a function of J/ ra( j (only for those galaxies shown in Figure|9p), are plotted 
in terms of 8/p as given in Equation |26| The ove rplot ted box corresponds to 
the physically acceptable ranges of 8 and p (see j)7.1> . 



global scale-length data, we plot the slopes of each fit in terms 
of 8 and p (as given in the exponent of Equationl26l> in Figure 
[TOl The slope of a given line in Figure [10] corresponds to a 
line in Figure [9J3 described by setting its slope equal to the ex- 
ponent in Equation|26] The physically acceptable ranges of 8 
and p are given by the overplotted rectangle. The steepest and 
shallowest gradients correspond to coordinates of the lower 
left and upper right corners of the rectangle, respectively. 

Only the fit to the structure scale-length data passes through 
the region of physically acceptable values for p and 8. In 
particular, nearly the entire range of p values, and values of 
8 > — I, appear to be acceptable. This suggests that diffu- 
sion occurring in star-forming regions is best characterized as 
having a significant contribution from soft scattering off of 
perturbations in the magnetic field. This scenario seems rea- 
sonable as the field near star- forming regions is k nown to be 
amplified and highly turbulent (Be ck et al. II19961) . 

The fit to the disk scale-lengths passes very near the plotted 
region of physically acceptable scaling values in Figure [10] 
Specifically, values of p = | (weak magnetic field structure) 
and 8 = — 1 (hard scattering) appear the most relevant. 

Finally, we find no physically acceptable combination of 8 
and p to reproduce the trend found between / g i b and (7 ra d. 
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This result, in addition to the order of magnitude estimate in 
37.1.11 suggests that systematic variations in ISM parameters 
alone cannot properly describe the mean distance traveled by 
a galaxy's CR electrons; rather, the age of the population must 
be the dominant parameter. 

7.2. Timescales of Star Formation Episodes 

Assuming that CR electron age is the dominant parameter 
which determines the measured best-fit scale-lengths, we can 
infer ages for the star formation episodes responsible for in- 
jecting the associated CR electron population into the ISM. 
We derive these ages from Equation[24] substituting / coo ] val- 
ues with the best-fit scale lengths. 

The mean equipartition strength of the total magnetic fiel d 
for typical spiral galaxies is found to be ^9 p.G dNiklas 1 1995b . 
In strong spiral arms, such as those of NGC 5194, the total 
field stre ngth is found to be s ignificantly larger at a value of 
~30 /jG dFletcher et al. 120 04). The field strength in more qui- 
escent regions of gala xies, like that near the Sun, have typica l 
values of ~6 /jG (e.g. lStrong. Moskalenko. & Reimer 1 12000). 
Let us assume that each of these characteristic magnetic field 
strengths correspond to the typical strengths for our global, 
structure, and disk components, respectively, and scale with 

1/2 

{/ d for the range in t/ ra( j found for the sample. 

The estimated ages for the global CR electron populations 
among the sample galaxies range from 0.14— 14. x 10 6 yr 
with a mean age of ^3.7 Myr. The ages for the CR electron 
populations associated with the disk and structure compo- 
nents range from 0.77 - 15. x 10 6 yr and 0.19 - 3.4 x 10 6 yr, 
and have mean values of ~5.3 and 1.2 Myr, respectively. Ap- 
propriately, we find the oldest ages to be larger than 10 7 yr; 
this is approaching the average cooling time of a 1.4 GHz 
emitting CR electron in a normal spiral galaxy. We also find 
that CR electrons residing in each galaxy's diffuse disk are 
>4.4 times older, on average, than the population of CR elec- 
trons associated with star-forming regions. 

Furthermore, the galaxy exhibiting the highest star forma- 
tion activity and shortest best-fit structure scale-length ap- 
pears to have a dominant CR electron population with an age 
of ~2x 10 5 yr. For this extreme case, such a young age 
suggests that the galaxy may still be in a phase of CR elec- 
tron generation within the shells of SNRs; assuming an am- 
bient particle density of 1 cm~ 3 and an explosion energy of 
10 51 erg, this age is ^5 times larger than the adiabatic life- 
time of SNRs, which by then should have expanded to have a 
diameter on the order of ~60 pc (e.g. ICondon II 1 992*1) . 

7.3. Morphologies of Star Formation and CR Escape 

In order to provide a physical interpretation for the location 
of galaxies in the best-fit scale-length — t/ ra d diagram, we de- 
scribe their general locations, designated by roman numerals, 
in Figure QT| So far in this discussion, we have not consid- 
ered the four lowest surface brightness galaxies (NGC 925, 
NGC 4725, NGC 4559, and NGC 3031) because their place- 
ment in Figure [9^ may be the result of a signal-to-noise ef- 
fect. Since it is difficult to determine whether this is in fact 
the case without having deeper imaging, we will now specu- 
late on a physical scenario which could explain their positions 
in Figure [9^. Similarly, the irregular galaxies in our sample 
(Holmb II, IC 2574, and NGC 4236) have also been excluded 
thus far in our discussion due to their lack of a disk compo- 
nent in the FIR and radio. The placement of these objects in 
Figure QT| and their behavior will now be put into the context 



of our phenomenological picture. 
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FIG. 1 1 . — The best-fit global scale-lengths for the entire 18 galaxy sample 
using isotropic kernels having an exponential profile and our free-free cor- 
rected radio maps. Overplotted on the scatter diagram are roman nume rals to 
identify the general placement of galaxies on this diagram (see j|7.3| for the 
full discussion). 



I. Peak of Most Recent Star Formation Event: In this re- 
gion of the diagram we assume a galaxy has just reached 
a peak in its surface brightness after a recent episode 
of enhanced star formation; the episode must have oc- 
curred within the last few Myr to ensure that SNRs are 
still relatively young and that the morphologies of the 
radio and infrared emission are similar because CR elec- 
trons have not had time to diffuse very far. While the 
nucle us of NGC 6946 has been classified as a mild star- 
burst dBall et al. 1 19851) . we note that our sample does not 
contain any systems for which the entire galaxy disk is 
experiencing a strong starburst. For systems experienc- 
ing such a mode of star formation, including the nuclear 
region of NGC 6946, high radiation fields and strong 
winds associated with the starburst could complicate this 
picture by rupturing magnetic field lines and increasing 
CR electron escape. 

II. Enhancement and Fading of Star Formation: After 
reaching a peak in global star formation activity, those 
CR electrons associated with the recent enhancement of 
star formation will dominate the CR population. As they 
spread through the galaxy they will begin to lose their 
energy to synchrotron and IC processes. In 37. 11 we have 
demonstrated that a galaxy's position along this trend 
cannot be due to variations in ISM parameters alone; 
rather, its star formation history is the dominant parame- 
ter. 

Alternatively, galaxies may move along this part of phase 
space as star formation activity becomes enhanced and 
then fades. Assuming a mean B of 10 p.G and £/ ra d = Ub, 
1 .4 GHz emitting CR electrons should lose their energy 
in ~1.9xl0 7 yr which, from Equation l24l corresponds 
to a propagation length of ^2.4 kpc. As this distance 
is consistent with the maximum values measured for our 
sample, we conclude that it takes a galaxy ~5xl0 yr 
to completely traverse this region if it stopped forming 
stars after a single episode and all Type II SNe occur on 
timescales <3xl0 7 yr. This is of course a lower limit, 
since galaxies have numerous star formation sites with a 
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FIG. 12. — We plot the global </fir ratio for the entire WSRT-SINGS sample as a function of FIR luminosity, radiation field energy density, and morphological 
type in panels a, b, and c, respectively. The horizontal line in each panel indicates the average q^iR ratio for the entire sample. Galaxies that are resolved at scales 
< 1 kpc, and therefore included in the two-component image smearing analysis, have been labeled. 



range of ages in their disks and star formation in spirals 
never completely turns off. 

III. FIR Heating by an Old Stellar Population ?: Galaxies in 
this region are characterized by low star formation ac- 
tivity and have best-fit scale-lengths which are shorter 
than expected given the trend found in region II. While 
we expect the best-fit scale-lengths to reach some fixed 
value for galaxies having constant SFRs over timescales 
comparable to, or longer than, the CR electron cooling 
timescale, and not increase indefinitely with decreasing 
infrared surface brightness, several possible explanations 
can be considered for the observed turnover. 

Let us assume this is not a signal-to-noise effect. One 
way to explain the shorter than expected scale-lengths 
is with an increase of CR electron escape. If this were 
the case we would then expect to find global FIR/radio 
ratios that are systematically higher than those for the 
more active star-forming galaxies. Alternatively, this 
observation could be explained if the less active star- 
forming galaxies have larger thermal radio fractions than 
the higher surface brightness galaxies. In this scenario, 
we would instead expect to find global FIR/radio ratios 
which are systematically lower than what is observed for 
the galaxies with higher radiation field energy densities. 

A non-linearity is in fact observed for the global FIR- 
radio correlation; galaxies having Lfir < 10 9 L© have 
radio luminosities which are lower than wh at their FIR 
luminosities would predict dCondonlll992l) . In Figure 
[P2h we plot the global ^fir ratios versus FIR luminos- 
ity for the entire WSRT-SINGS sample. The known 
non-linearity is clearly observed. In Figure [12b we in- 
stead plot ^fir against t/ ra( i; here we do not see any ev- 
idence for the lowest surface brightness galaxies to ex- 
hibit systematically larger or lower FIR/radio ratios. It 
then seems that the above explanations (i.e. increased 
CR electron escape or larger thermal radio fractions) are 
not able to reconcile the short scale-lengths found for 
objects in region III. 

We also note that the trend found between the best-fit 
structure scale-lengths and radiation field energy den- 
sity does not exhibit a turnover from regions II to III 
as seen for best-fit disk and global scale-lengths. This 
suggests that the propagation of CR electrons associated 
with star-forming complexes is largely similar despite 



the differences in the galaxy's global star-formation ac- 
tivity. Then, to simultaneously explain this result from 
our two-component analysis along with the depressed 
best-fit global scale-lengths, we may consider a situa- 
tion in which older stars contribute significantly to the 
FIR heating. This will make the FIR emission appear 
more diffuse and require less smoothing to best match 
the distribution of radio emission, as observed. 

IV. Unpopulated Part of Phase Space: We find no galaxies 
occupying this region of phase space in which relatively 
short global scale-lengths would be measured for a mod- 
erate values of f/ ra( j. One way for a galaxy to populate 
this part of the diagram is if, after a long period of qui- 
escence, it begins to form stars at a moderate rate. The 
galaxy would then pass through this region very quickly 
(<10 7 yr) before shifting into region II. The lack of such 
galaxies suggests that, at least in the local Universe, star 
formation in spirals does not completely cease for long 
periods of time. 

V. CR Electron Escape: As stated in 36.11 the irregulars 
behave markedly different than the sample spirals. We 
do not believe this discrepant behavior to be the result 
of signal-to-noise effects. However, we do note that 
the radiation field energy densities of these galaxies, as 
derived from their infrared surface brightnesses alone, 
may be significantly underestimated due to the low dust 
content and more transparent media in irregular galax- 
ies. While shifting these objects to larger values of <y raa - 
in Figure QT| they would still clearly remain outliers 
compared to the rest of the sample. The morphologies 
and ISM of these galaxies do not seem consistent with 
keeping their CR electrons bound outside of their ini- 
tial clouds around SNRs. After the CR electrons leave 
this cloud, the lack of a dense ISM and magnetic field 
to keep them trapped through multiple scatterings off 
of magneto-hydrodynamic (MHD) waves allows them to 
easily escape the system and enter intergalactic space. 

Evidence for this picture comes from the extremely short 
best-fit scale-lengths measured for these galaxies. The 
inability to retain CR electrons once they leave dis- 
crete star-forming regions should lead to deviations in 
the FIR-radio correlation; ratios should be significantly 
higher than the nominal value which is in fact observed. 
In Figure [12b we plot the global FIR/radio ratio as 
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a function of morphological type. All three irregular 
galaxies in our sample clearly exhibit FIR/radio ratios 
which are higher than the average value for the entire 
sample; the average qpi& for the irregulars is ^0.34 dex 
larger than the mean for the entire sample. Further- 
more, detailed multi-wavelength studies of the SINGS 
dwarf galaxies have uncovered these deviations within 
individual systems as well. Strong variations (~1 dex) 
have been found in the TIR/radio ratio of IC 2574 
dCannon et al. 1 120051) while the radio continuum emis- 
sion in the dwarf starburst NGC 1705 has been found 
to be an order of magnitude less than w hat one would 
expec t from the nominal TIR/radio ratio (C annon et al. I 
20061). A more detailed, systematic study of the FIR and 
radio properties for dwarf galaxies is beyond the scope 
of this paper. 

8. REMAINING ISSUES AND FUTURE PROSPECTS 

At the present, we appear to be reaching the limitations 
on what we can learn about the evolution of CR electrons 
in nearby galaxies with current data. Using more reli- 
able SFR tracers (e. g. 24 ^m + Ha; ICalzetti et al~l 120071; 
iKennicutt et al. Il2007l) to compare with the radio should im- 
prove our ability to understand the effects discussed in 37.31 
Modeling is also needed and should lead to quantitative esti- 
mates for the SFR variations from the placement of a galaxy 
in Figure [TT] this can be a powerful tool for constraining the 
recent star formation histories of these galaxies. 

Due to our small sample size, especially for the lowest sur- 
face brightness galaxies, a number of outstanding issues re- 
main. For one, it is currently difficult to determine whether 
the turnover seen in Figure [TT] is in fact physically driven or 
simply due to a signal-to-noise effect; deeper integrations are 
clearly needed. At the very high surface brightness end of the 
spectrum there is even less information to gauge how well our 
proposed phenomenology applies since the present sample 
lacks any starbursting luminous infrared galaxies (LIRGs); 
such galaxies are typically too distant to resolve on < 1 kpc 
scales with current FIR telescopes. We speculate that such 
galaxies will behave similarly to the high surface brightness 
galaxies in our sample, although strong galactic winds aris- 
ing from the central starburst could complicate the picture by 
transporting CR electrons directly out of the disk, or perturb- 
ing severely the magnetic field structure. The fact that such 
galaxies obey the FIR-radio correlation suggests that the mag- 
netic field must be amplified for synchrotron cooling to retain 
the same fraction of CR electron energy losses relative to es- 
cape and other competing processes. It is only with additional 
observations that these outstanding issues can be addressed. 

Applying the smoothing phenomenology to such galaxies 
will add important new insight on the role of magnetic fields 
and CRs in the process of star formation and galaxy evo- 
lution as a whole; starbursting ultra-luminous galaxies con- 
tribute significantly to the luminosity density of the Universe 
during the most critical epochs (z > 1) of galaxy evolution. 
More than 90% of UV and visible radiation from starbursting 
ultra-luminous LIRGs (ULIRGs) is absorbed by dust leaving 
FIR and radio observations as the most promising way to gain 
insight on the star formation processes within such systems 
and better understand the physical connection between the 
gaseous and relativistic phases of the ISM. Such a connection 
includes that of feedback processes; the production and expul- 
sion of CRs in galaxies ma y work to cap their luminosity and 
star formation intensity ( Socrates, Davis, & Ramirez-Ruiz 



120071) . The outlook for such additional observations appears 
rather bright; with Herschel and the EVLA on the horizon, 
increasing the sample size and improving the sensitivity and 
resolution of the observations allowing the inclusion of more 
distant galaxies should become achievable within a few years. 

9. SUMMARY AND CONCLUSIONS 

Using a two-component image-smearing analysis, we have 
separated the signatures of CR electron diffusion at spatial 
scales corresponding to star-forming structures (<1 kpc) and 
galaxy disks (>1 kpc) within 18 galaxies observed as part of 
SINGS and WSRT-SINGS. Our results and conclusions can 
be summarized as follows: 

1 . We confirm and extend earlier results of M06a,b. Em- 
pirically, the dispersion in the FIR-radio correlation 
within galaxies is most reduced by an image-smearing 
model; this improvement is significantly better than 
what can be achieved by fitting correlations and remov- 
ing linear trends. 

2. The best-fit global scale-lengths decrease as a function 
of increasing star formation activity as measured by the 
infrared surface brightness of a galaxy. Our interpreta- 
tion is that a galaxy's CR electrons are closer to their 
place of origin within galaxies having intense star for- 
mation activity. 

3. The trend of decreasing best-fit global scale-length with 
increasing radiation field energy density is due to higher 
surface brightness galaxies having undergone a recent 
enhancement of star formation rather than variations in 
other ISM parameters. For sufficiently large enhance- 
ments, these galaxies are observed within ~40 8 yr of 
the onset of the most recent star formation episode. 

4. Unlike spirals, irregular galaxies lack any well defined 
diffuse disk component at either 70 /jm, or especially 
at 22 cm. Presumably, the CR electrons escape these 
galaxies soon after leaving their parent star-forming re- 
gions due the absence of a dense ISM which would 
keep large-scale interstellar magnetic field locked into 
place. This conclusion helps to explain why these 
galaxies have global FIR/radio ratios systematically 
greater than canonical values. 

5. As infrared surface brightness increases, the character- 
istic diffusion scale-length of a galaxy's CR electron 
population begins to transition at log f/ ra( i < —12.5, or 
log Esfr < —2.3, from being biased by CR electrons 
making up its diffuse disk to being biased by those re- 
cently injected near star-forming structures. From this 
we conclude that a galaxy's CR electron population 
transitions from being dominated by old CR electrons 
to being dominated by young CR electrons as a function 
of star formation intensity. 

6. The two-component analysis works better than smear- 
ing with a single smoothing kernel for spiral galaxies 
of type Sb or later which have high amounts of ongo- 
ing star formation activity (i.e. ^40% of the sample). 
This result suggests that star formation must be intense 
and highly structured for the two-component analysis 
of these data to differentiate properly between the dif- 
ferent CR electron populations. 
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